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FOREWORD 


This report describes the results of pogo stability investigations under- 
taken at The Aerospace Corporation during the one-year period beginning 
October 1975. These investigations were a continuation of the studies reported 
in Ref. 4 under NASA contract NAS9-14142; the NASA(JSC) Technical Monitor 
was Dr. H.H. Doiron. The authors would like to acknowledge the work of 
R.E. Orth in programming support for the studies, the contribution to the 
transfer function error analysis by Dr. J.S. Bendat, and the constructive 
comments of Dr. Doiron during the course of the investigation. 
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ABSTR.^Cr 


Another phase of the study of pogo suppression during Space Shuttle 
Main Engine (SSME) operation was conducted as a follow-on to the study 
reported in Ref, 4. The latest available structural modal data (from model 
5, 3A) were incorporated and stability studies with the planar multiengine 
model were conducted for the lox system at the liftoff time of flight. The 
sensitivit ' of the system stability was determined selectively for dispersed 
values of certain cri^^ical system parameters on the basis of previous studies 
and recent experimental findings. The dispersions involved were (1) toler- 
ance on structural natural frequencies^ (2) increased downcomer axial motion 
in the 2, 54-Hz structural mode along with high pump gains, (3) decrease' * 
resistance of the low pressure oxidizer turbopump (LPOTP), and (4) a large 
upward shift of the frequency of the upper SSME structural mode to about 
30 Hz. In all cases the presence of the baseline accumulators virtually 
eliminated the destabilizations occurring in the absence of accumulators. 

Studies of the stability consequence of reduced accumulator compliance 

at liftoff yielded a minimum recommended compliance equivalent to that 

3 3 

provided by a gox volume of 0.007 m (0. 25 ft ). A limited study of the con- 
sequence of accumulator failure showed that loss of function of the upper 
SSME accumulator led to instability under toleranced conditions in the shifted 
upper SSME structural mode. 

Although not investigated in this study, it should be noted for complete- 
ness that previous investigations have identified a potential for instability in 
the separation mode introduced by the accumulator (fluid oscillation between 
accumulator and the high pressure oxidizer turbopump (HPOTP). These 
earlier studies showed that accunmlator resistance can be highly effective in 
eliminating such an instability. 
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Study of the refinement of the analytical model of the LPOTP was 
initiated by investigation of curve fitted transfer functions and by the formu- 
lation of a first-cut physical model, all based upon recent data from dynamic 
testing of a quarter-scale model impeller. Although both approaches showed 
promise, certain questions remain to be resolved before an implementation 
into the system stability model can be made with confidence. Further studies 
are planned. 

A preliminary sensitivity analysis for pogo stability during operat r. of 
the orbital maneuvering subsystem (OMS; was performed. Included were the 
development of OMS pod structural and propulsion system models. The 
results of the study revealed a strong destabilizing interaction in a high-gain 
engine axial mode, as well as a potential for other destabilizing interactions. 
An initial evaluation indicated the possible benefit of small accumulators. 

Investigations were made of various sources of error in the experi- 
mental measurement of frequency response functions as a result of system 
noise. An error estimate was developed to account for r'lndom errors due to 
the noise, with application to both fixed and swept frequency sinusoidal test- 
ing. Application to the planning of a test was exemplified. The results are 
meaningful for the planning and* interpretation of dynamic response from 
propulsion tests. 
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1. POGO SUPPRESSION FOR MAIN PROPULSION 
SUBSYSTEM OPERATION 


This section of the rt'port describes the continuation of the stability 
studies undertaken by The Aerospace Corporation in support of pogo sup- 
pression during operation of the main propulsion subsystem (MPS). Earlier 
studies presented initial analyses of the Shuttle pogo susceptibility (Ref. 1) 
and the use of single equivalent engine stability models to investigate the 
effectiveness of passive and active suppressors (Refs. Z and 3, respectively). 
These studies were followed by the development of a multiengine pitch-plane 
stability model (Ref. 4) under the initial phase of the current contract with 
the NASA Johnson Space Center (JSC). 

Included in this section are an outline of the analytical model and 
numerical procedures employ<;d to determine* the system stability, a discus- 
sion of the numerical input data for the structure and the propulsion system, 
and the results of the updated stability analyses. 

1. 1 ANALYTICAL MODEL AND NUMERICAL PROCEDURES 

The launch configural ion of the Space Shuttle, shown in Figure 1, 
comprises the orbiter vehicle, and extcirnal tank, and two solid rocket 
boosters. The external tank contains a forward liquid oxygen tank and an 
aft liquid hydrogen tank for supply of the three SSME*s located on the orbiter 
vehicle. The relative location of the lox tank and the SSME's (see schematic 
in Figure 1) necessitates an extensive length of feedline. The total length is 
about 40 m (130 ft), which introduces several propulsion system modes in 
the frequency range of interest. 

The multiengine pitch-plane model developed in Ref. 4 was employed 
as the basis for the present system stability analysis. The model comprises 
the vehicle structure, the SSME's and the elements of the MPS associated 
with the lox circuit. The hydrogen system does not appear in the model 
since early assessment of the pogo problem indicated that the hydrogen system 
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Figure 1, Schematic of Space Shuttle and Lox Circuit 



Figure 2. Schematic of Multiengine Model 
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was a less likely contributor to potential instability. The geometric arrange- 
ment of the model is shown in Figure 2 and the associated system equations 
are presenteu in Appendix A. Observe that only one lower engine is illus- 
trated in the figure since the motion of both lower engines is identical for this 
pitch-plane model. 

To assure a good description of the higher organ-pipe modes of the 
lox feedline, the flow in the long feedline segment that is located along the 
external tank (the element between points and in Figure 2) was repre- 
sented by the exact solution for one - dimensiona 1 continuous compressible 
flow (Ref, 1), The fluid flow in the remaining feedline segments and in the 
inter-pump ducts was describ«;d by a combination of incompressible flow seg- 
ments and local compliances to represent liquid compressibility corrected 
for pipe radial elasticity. The compliances were located at feedline corners 
and at the inlets and outlet® of the low and high pressure oxidizer turbopumps. 
The resulting flow description provided an adequate representation of distri- 
buted compressibility effects in these line elements for frequencies up to 
30 Hz. Comparison of the lox system modal frequencies predicted in this 
manner with those determined from a fully compressible model showed ac- 
ceptable agreement for this range. 

The dynamic behavior of the LPOTP's and the HPOTP's was described 
by the following equations relating the downstream pressure P^, and flow Q^, 
to the upstream pressure P^, and flow Q^. 


Pj = (m + l)P^ - (Ls + R)Q^ 


Q 



sC P 
P ^ 


( 1 ) 


where L, R, and m+1 denote frequency independent inertance, resistance, 
compliance, and pump gain, respectively; s is the Laplace variable used to 
denote the complex frequency. This "elementary" pump model has been used 
in the pogo stability analysis of Titan, Thor, Delta, Atlas and Saturn boost 
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vehicles (Ref* 5). Until the tests performed on the LPOTP model impeller, 
by Caltech in 197^^ (Ref. 6), there w.is no experimental basis for a more 
refined representation of the dynamic behavior of the cavitating turbopump. 

The accumulators used for pogo suppression were represented by 
an inlet liquid inertance and a gas < ompliance. The program also includes a 
capability for the input of accumulator resistance; however, this capability 
was not employed in the present an,tlyses. For these studies, the accumula- 
tor is located two line diameters u]>stream oi the inlet flange of each HPOTP. 
This separation is significant in that it both influences the effectiveness of 
the accumulator and introduces a n<^w propulsion system mode of vibration, 
termed the "separation” mode, that provides an additional possibility for 

instability (Refs. 4 and 7). This, particular separation represents a worst 

❖ 

case interpretation of the Shuttle system requirements. 

The computer mechanization of ihe analytical model of the system 
allowed for representation of the motion of the Shuttle vehicle by a series of 
up to ten pitch-plane structural modes. Tank dynamic outflow effects are 
included in the generalized force terms on the basis that the structural modes 
used in the analysis were developed with closed-bottom tanks (Ref. 1). The 
eigenvalues and eigenvectors of the coupled structural/propulsion system 
were developed from the condition that the determinant from the equations 
for free vibration of the coupled system vanishes. Transcendental functions 
appear within the determinant because of the exact distributed representation 
of the downcomer flow and the pro\ision for arbitrary phase shifting of the 
engine thrust forces. The solutions of the characteristic equation were 
determined with the use of an iterative root-finding subroutine {Ref. 8) 
that used the input structural mode frequency and previously calculated 
propulsion system frequencies as initial guesses. 


"I' 

The system requirements called for the accumulator to be located no more 
than 0. 33 m (1^ in. ) from the pump inlet flange (approximately two line 
diameters); the actual location is 0. 17 m (7 in. ) from the inlet. 
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NUMERICAL INPUT DATA 


The structural modal data xised in the analysis consisted of frequen- 
cies and amplituu 3 of a selected set of pitch plane modes developed from 
Shuttle model 5.3A and provided by Roc kwell International/Space Division. 

The earlier model 4. 3 data were employed in the previous study (Ref. 4). The 
modes selected by Rockwell were based upon the magnitude of the modal 
amplitudes for those elements of the system that were of concern in the pogo 
problem (i. e.^ engine gimbals, feedlinc* corners, etc. ) for frequencies up to 
30 Hz. Although such modal sets were provided for seven flight conditions, 
this analysis dealt only with the liftoff e vent. In the assignment of modal 
amplitudes, the corners immediately upstream of the LPOTP and HPOTP 
positions (points^), and in Figure Z) were assigned the same 

motion as the associated pump. The motion of the manifold and the longitudi- 
nal displacement of the corner immediately ahead of the manifold (points 
and^^in Figure 2) were assumed to be zero; this assumption was necessi- 
tated by the fact that the strudural definition was not sufficiently adequate to 
define these modal motions. The assumption was not considered to be critical 
since the length of the feedline segment between points(^and^)was small 
compared to the length of the downcc>mer. This relative insensitivity had also 
been confirmed by results in Ref. 4. Examples of the structural mode shapes 
employed in the analysis are given in Appendix B, where the data for the 
second and fifty-fourth modes at liftoff are provided. These modes proved 
to be of particular interest from the stability analyses. To account for 
damping in the vehicle structure, a critical viscous damping ratio of 0.01 was 
assigned to each structural m(*de. 


The cross-sectional area, length, and radial elasticity of the various 
lines were based upon available design data. The resistance and inertance of 
the lines, pumps, engine, and thrust chamber were developed from the basic 
data given in the SSME Engine Dynamic Model (Ref. 9). The values of these 
parameters are provided in Appendix B, together with the local compliance 
values employed in the feedline description. Values for the pump cavitation 



compliances and their uncertainty ranges as a function of net positive suction 
head (NPSH) were based upon the i^mpirical studies of Ref* 10 and applied to 
the Shuttle pvux.^^s in Ref. 1. The pump gain, m + 1, was related to the com- 
pliance based upon an empirical relationship developed from unpublished 
Titan and Delta vehicle studies* 

m+1 = 1 + 3000 pg Cp/A^ (2) 

where p is the liquid density, g is the gravitational constant, and A. is the 
flow area at the pump inlet. The variations with flight time of the estimated 
nominal values of the compliance and gain parameters are illustrated in 
Figures 3 and 4* The source of the operating pressures versus time is the 
Space Shuttle Pogo Data Manual (Ref. 11). The nominal, maximum and mini- 
mum values of compliance, and the corresponding pump gains for the liftoff 
event are provided in Table 1. 

1.3 PROPULSION SYSTEM MODES 

Knowledge of the propulsion system mode ^ not only provides initial 

guesses for the numerical determination of the eigenvalues of the coupled 

system, but also facilitates increased understanding of the coupling processes 

revealed by the stability analysis. Vibration modes of the propulsion system 

3 3 

were determined with and without the baseline 0*017 m (0,6 ft ) gaseous 
oxygen (gox) accumulators. The calculations were undertaken using the 
Aerospace nominal values for the pump parameters. The frequencies calcula- 
ted for the liftoff condition are shown in Figure 5, together with the frequencies 
of some of the structural modes of interest; the numbers of some of the struc- 
tural modes are also provided in the figure. Examination of the associated 
mode shapes of the propulsion system indicated that certain distinct Lypes of 
modes could be identified. 

One stich type was an "inter-engine" mode wherein the predominant 
modal perturbations took place within the engines and the feedline run between 
the engines. Two of the propulsion system modes calculated in the frequency 
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Table 1. Values for Pump Gain and x-nimp Cavitation 
Compliance at Liftoff 


Parameter 

Value 

LPOTP Gain, m^ + 1 

1.48/2.23/1. 10* 

HPOTP Gain, m 2 + 1 

1.63/3.64/1. 24 

IjPOTP Compliance, Cpj 


or? 0 

10 ' m /(MN/m ) 

1. 0/2. 4/0. 18 

-2 2 
(10 in.'^) 

(1.6/4. 0/0.3) 

HPOTP Compliance. Cp 2 


10"^ m^/(MN/m^) 

0. 36/1. 5/0. 14 

-2 2 
(10 in/) 

1 

(0.60/2.3/0. 23) 


Nominal/maximum/minimum values provided for all 
entries 

Compliance in SI units is a volume change per urjit 
pressure change; in English units, it is the weight of 
liquid displaced per unit pressure change. 
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range of interest had been identified as being cf this type and the frequencies 
are designated with an l/E notation in Figure 5; the distribution of fluid 
amplitudes associated with the first of these modes is illustrated in Figure 6, 
As noted in the figure, the fluid flow is identical in the two lower branches 
(the figure only shows one of the lower branches). In the ^ase of a system 
without accumulators it will be note<l from the results prt ented in Figure 5 
that the addition of the baseline accumulator has a dramatic effect on the 
frequency of the tirst inter-engine mode, the frequency being reduced from 
about 14 Hz to around 3. 5 Hz by the addition of the accumulator. The fre- 
quency of the second inter-engine mode is only slightly affected by the ac- 
cun v^lator. 


Another mc.ie was identified as an "inter-pump** mode and the asso- 
ciated frequency is designated I/P in Figure 5. This particular mode involves 
a considerable amount of fluid motion between the compliances provided by 
the UPOTP and HPOTP cavitation. Although the overall mode shape did not 
exhibit a pure inter -pump character, the mode was identified in this manner 
since the frequency was estimated reasonably well by the following formula 
for pure interpump action: 


f. 

1 


1 

ITT 


1/C ,+1/C 



L, + L , 

1 pi 



1/2 


(3) 


where L., denote the inertance of the inter-pump line and of the LPOTP, 
resp>ectively; and where cavitation compliance values at 

the inlet to the LPOTP and HPOTP, respectively. The incorporation of the 
baseline accumulator reduced the frequency cf this mode from about 19 Hz 
to around 11 Hz. The remaining propulsion system modes showm in Figure 5 
^have shapes associated vdth acoustic modes of the feedline from the tank. 

The frequencies of these acoustic modes lie within the frequency bands de- 
fined by the open-open c»nd open-closed modes of the feedline, and there is 
relatively little participation within the engines. This character is illustrated 
in Figure 7, where the frequencies and associated damping ratios for the 
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Figure 6. Amplitudes in the First Inter-Engine Fluid Mode 

at Liftoff 
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propulsion system modes are given for the cases of a system with;#ut accu* 
mulators and a system with the baseline accumulators. The figure shows 
that the higher fcedline modes are lightly damped and that the inter-pump 
and inter-engine modes are relatively highly damped. Figure 7 shows that 
the frequency reductions in the feedline-type modes, due to the presence of 
the baseline accumulators, becomes irsignificant above the lowest two modes. 
This behavior occurs since the presence of the LPOTP compliance is serving 
to provide a nearly open boundary condition for the higher feedline-type 
modes. 

Finally, it should be noted that the frequency of the separation mode 
(i.e*, the propulsion system mode associated with the fluid between the 
acciimulator and the HPOTP inlet cavitation bubble) lies above the frequency 
range of interest for the Aerospace nominal values of the pump parameters 
(namely, 38 Hz). The effect of this mode on the system stability was not 
investigated in the present study. However, it should be remarked that 
previous investigations (Refs. 4 and 7) have shown that this fluid mode rep- 
resents another potential source of instability. These studies indicated that 
accumulator resistance is highly effective in eliminating such instabilities. 

1.4 STABILITY ANALYSIS RESULTS 

The stability analyses focused upon the liftoff event. A total of 
thirteen structural modes was analyzed to determine the stability situation at 
this flight condition. The analyses were conducted with one structural mode 
present at a time, and the input structural mode frequency was dispersed 
^15 percent about its nominal values to account for uncercainties in the 
structural modeling. This variation was judged to be a reasonable estimate 
to cover the worst case conditions in terms of the proximity of the structural 
and propulsion system resonances (see Ref. 4). The damping ratio used for 
all structural modes was 0.01. The analyses were generally conducted using 
the Aerospace nominal values for the parameters appearing in the elementary 
pump models for the LPOTP and HPOTP; the exceptional case was the use of 
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the Rockwell nominal values for the pump cavitation compliance in V»o 

the sensitivity studies that were conducted with a **shifted” structural mode 

frequency (see Section 1.4.6). The system without suppression and the 

3 3 

system with the baseline 0.017 m (0.6 ft ) gaseous oxygen (gox) accumu4.c. j 
were first treated. These nominal case analyses were then followed by a 
series of sensitivity studies in which the effects of downcomer motion, 
pump gain, accumulator volume, LPOTP resirtance, structural mode fre- 
quency shift, and accumulator failure were examined. 

1.4.1 Nominal System; No Suppression 

In the study of Ref. 4 the liftoff stability analys ;s undertaken with 
the multi-engine model and the model 4,3 modal data indicated that the 
nominal system was free of pogo instability. The present analysis, under- 
taken with the model 5. 3A set of modal data, indicated that the stability of 
the system had been degraded relative to the previous results. Specifically, 
it was found that 

a. The L2 mode at 2. 54 Hz gave rise to a minimum damping 
ratio of 0.0053, whereas the previous minimum damping 
ratio was 0. 0087. 

b. An instability was now predicted in a new local upper SSME 
structural mode at 17.8 Hz ^L54 mode). 

The stability behavior in the L2 mode is illustrated in Figure 8 
where the results obtained for both the 4.3 and 5. 3A modal sets are pre- 
sented. The destabi zing influence in this mode results from coupling be- 
tween the structural mode and the fundamental mode of the propulsion system 
(see Figure 5), The source of the degradation in system stability is the 
change in structural mode shape. Specifically, the axial modal amplitude of 
the downcomer is larger in the later modal data. Increase in modal ampli- 
tude at this location was previously shown to reduce the stability of the L2 
mode (Ref. 4). 

The L54 modal amplitude distribution is illustrated in Figure 9. 

The sensitivity of stability in this mode to structural mode frequency is 
f shown in Figure 10, where it is seen that a minimum damping ratio of -0. 002 
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Note; Upper engine motion predominates. 
Arrows not to scale. 



Figure 9* Modal Amplitudes in L54 Mode 


Figure 10. System Damping without 
Accumulators: L54 Mode 





is predicted at a structural frequency of about 14 Hz, Examination of the 
eigensolution for this case reveals that the instability results from coupling 
of the structural mode with the first inter-engine mode of the propulsion 
system (Figure 5), Significant coupling takes place because of the similar 
shape of the involved structural and propulsion system modes (compare 
Figures 6 and 9). 

1.4.2 Nominal System; Baseline Accumulator 

3 3 

The effect of the baseline 0. 017 m (0. 6 ft ) gox accumulator on the 
frequencies of the propulsion system modes was discussed in Section 1.3, In 
particular, it was seen from the spectra presented in Figure 5 that the fre- 
quencies of the fundamental mode and of the first inter-engine mode were 
decreased by some 40 and 75 percent, respectively. is not surprising then 
that the stability in the L2 mode is enhanced (see Figure 11) and the instability 
in the L54 mode eliminated (see Figure 12) by the presence of the accumula- 
tor. One question which does arise, however, is whether the inter-engine 
mode, which is now the second mode of the propulsion system (see Figure 5), 
significantly couples with the lower structural modes. This question was 
examined for the L2 mode and it was found that no significant destabilizing 
coupling took place. The reason for the weak coupling is that the inter- 
engine mode has fluid motion in the upper engine which is out of phase with 
that in the lower engines, whereas the structural motion of the engines are 
all in phase- consequently, a cancelling effect takes place. Examination of 
the modal character of other low frequency structural modes indicated that 
significant coupling with the inter-engine mode would be unlikely. 

1.4.3 Downcomer Motion /Pump Gain Study 

The stability of the L2 mode had previously been shown to be sensi- 
tive to both the modal amplitude of the downcomer and pump gain (Ref. 4). 

This sensitivity was checked with the present modal data with and without 
accumulators. Analyses were run with a range of downcomer modal ampli- 
tudes for unit pump gain, the Aerospace nominal, and the Aerospace maxi- 
mum pump gains. The resulting minimum damping ratios, calculated over 
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Figure 11. Effect of Accumulators on System 
Damping: L2 Mode 


Figure 12. Effect of Accumulators on 
System Damping: LS4 Mode 


J 


( 


I 


> 

i 


the specified frequency range (2. 54 Hz i: 15 percent), are shown in Figure 13 
for the case of no accumulators. The damping levels are shown as a function 
of the ratio of the downcomer modal displacement to the average axial 

modal displacement of the SSME gimbals (f> (g). The results indicate that 

pvimp gain increase remains as a significant destabilizing influence, the effect 
tending to be more pronounced at the higher levels of downcomer motion. 

The character of the curves in Figure 13 reflects a balance between the de- 
stabilizing effect of engine thrust [increases with pump gains and with 
0 (g)]and the stabilizing effect of the axial pressure-area force at the aft 
corner of the downcomer [increases with 

creases, the stabilizing effect eventually dominates. The result is the appear- 
ance of a minimum, as seen in the curves for unit and nominal pump gains. 
Such a minimum would also occur in the curve for maximum pump gains if 

the plot extended to large enough values of 0 

X / xa 

The sensitivity of the L2 mode stability to pump gains is further 
illustrated by the results show.i in Figure 14. There the minimum damping 
ratios, for an arbitrarily selected worst case value of downcomer modal 
amplitude ratio of 0. 75, are presented as a function of the piomp gain product 
(mj + l)(m 2 + l). Results are given for the system with and without accumula- 
tors. The extreme sensitivity of the damping level to pump gain in the 
absence of suppression is clearly shown; it is seen that a 50 percent increase 
of the pump gain product from unity to a value of 1. 5 produces a 40 percent 
drop in the predicted damping level. The sensitivity is considerably reduced 
by the presence of the baseline accumulators. In this instance the damping 
ratio has only dropped to a value of about 0.007 for the maximum value of 
pump gain product (8. 1) that was treated; in the corresponding case with no 
accumulators, the system has become unstable with a damping ratio of about 
-0. 006. 

1,4,4 Accumulator Volume 

The sensitivity of the system stability to accumulator volume was 
also examined to establish a minimum required accumulator precharge at 
liftoff. This information was n supjjort of studies of alternative charging 
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Figure 13. Minimum System Damping versus 
Downcomer Motion for Different 
Pump Gains; L2 Mode 


MINIMUM DAMPING RATIO, 



Figure 14, Minimum Sy-»tem Damping versus 
Pump Gain Product; L2 Mode 



procedures for the accumulators. The study was directed to the L2 mode 
since previous analyses had indicated that volume sensitivity was more pro- 
nounced in the lower frequency modes (Ref. 4). The results of a series of 
analyses for a downcomer amplitude ratio of 0. 75 are presented in Figure 15; 
the minimum damping ratio is plotted against accumulator volume fox unit and 
nominal pump gains. The latter case shows a much greater sensitivity to 
accumulator volume and is therefore used as a basis for the recommended 

minimum accumulator volume at liftoff. It is seen that stability is strongly 

3 3 

enhanced as the accumulator volume is increased to about 0.014 m (0.5 ft ); 
a much lower rate of stability enhancement occurs for further increases in 
accumulator volume. 

To establish a minimum accumulator compliance that* should be 

available prior to liftoff, it was assumed that the accumulator was required 

to limit system stability loss to less than half of the structural damping. This 

was based upon the 6 dB damping gain margin recommended for nominal 

system parameters in Ref. 5. Thus, a system damping ratio of greater than 

0,005 is required for an assumed structural damping ratio of 0,01. The 

results in Figure 15 (for nominal pump gains) show that an accumulator com- 

3 3 

pliance equivalent to that provided by a gox volume of 0.007 m (0.25 ft ) is 
the minimum required prior to liftoff. 

1,4.5 LPOTP Resistance 

A study was directed at the sensitivity of the L54 mode stability to 
variation of the resistance of the LPOTP. Previous analysis had indicated 
that the L2 mode stability was insensitive to such variation. The incentives 
for this particular study were twofold. First, on a physical basis it was 
anticipated that the damping in the inter-engine mode of the propulsion sys- 
tem would be sensitive to LPOTP resistance. Second, the experimental 
results developed by Caltech for the LPOTP impeller (discussed in Section 
2, 1. 1) suggested that the resistive impedance decreased dramatically at high 
frequencies of oscillation. The results of an analysis undertaken for the 
case of zero LPOTP resistance revealed that loss of this resistance signifi- 
cantly degraded the stability of the system. Results are shown in Figure lb 
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Figure 15. Effect of Accumulator Volume on 
Minimum System Damping: L2 
Mode 
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Figure 16. Effect of LPOTP Resistance on 
System Damping; L54 Mode 



for Ihe system withou' accumulators for zero LPOTP resistance and for the 

5 ?. 

nominal resistance of 8. 8 MNs/m (0. sec/in. It is seen that the lowest 
predicted damping level has dropped from -0, , JZ5 to -0, 0095, Fortunately, 
the frequency shift efiect provided by the baseline accumulator is sufficient 
to desensitize the stability to LPOTl^ resistance. This is illustrated in 
Figure 16 where the results obtained with the baseline accumulator and zero 
pump resistance are also presented. These results are quite close to those 
obtained for the nominal LPOl'P resistance. 

1.4,6 L54 Mode Frequency Shift 

Subsequent to the performance of the stability analyses with the 
model 5. 3A modal data, some preliminary results from the horizontal ground 
sribration test (HGVT) that had been performed on Orbiter Vehicle 101 became 
available (Ref, 12). These results indicated that a local upper-engine mode 
occurred at about 30 Hz, rather than at 18 Hz as predicted in the modal 
analysis. Lacking mode shape data for this 30-Hz mode, it was decided to 
examine the effect of the new frequency by simply shifting the frequency of 
the L54 mode (a local uppex - engine mode) to the 25-32 Hz range. From the 
frequency spectra shown in Figure 5, it is seen that such a shift could bring 
this structural mode into proximity with the second inter-engine mode of the 
propulsion system with or without the accumulators present. For conserva- 
tism, the LPOTP resistance was set to zero. The results of analyses under- 
taken for the system both without accumulators and with the baseline accumu- 
lators are shown in Figure 17, In the absence of accumulators it is seen that 
there are no strong destabilizing tendencies. In fact, the opposite is true in 
the 27-29 Hz range where the system damping is significantly enhanced from 
the prescribed 0.01 structural damping ratio Examination of the generalized 
force components indicates that this enhancement is primarily due to the 
engine thrust forces. When the accumulators are introduced the thrust 
forces, and hence their stabilizing influence, are reduced. This effect is 
seen in the relationship of the stability curves for the case with and without 
accumulators in Figure 17. A similar effect was observed in previous 
studies (Ref. 4, Section 2. 3.2,2, Figure 27). 
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Figure 17. Stability of L54 Mode Shifted to 
Higher Frequency 
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Figu**e 18. Stability o^ L54 Mode Shifted to 
Higher Frequency: Rockwell 
Nomi.ial Pump Compliances 



It was also of interest to see what happened if the fi rst inter-engine 
mode frequency was shifted to this higher range. Such a shift could be 
achieved by replacing the Aerospace nominal pump cavitation compliance 

values of 0. 0012 m^/MN (0. 02 in.^) for the LPOTP and 0. 00027 m^/MN 

2 

(0.0045 in. ) for the HPOTP with the corresponding R< jkwell nominal values 
of C. 39 10"^ m^/MN (0, 0067 in.^) for the LPOTP and 0. 69 10“^ m^/MN 

7 

(0.001 in. ) for the HPOTP. The pump gains were not changed. St. .lity 
analyses undertaken for these conditions and no accumulators revealed that 
significant destabilizing coupling occurred between the L54 mode and the first 
inter-engine mode at these higher frequencies. These results are shown in 
Figure 18 where it is seen that the system is only marginally stable. It was 
this same type of interaction which had induced instability with the unshifted 
L54 mode and the Aerospace pump compliances (see Figure 10). The results 
obtained when the baseline accumulators are added to the system are also 
shown i’^ Figure 18 for the cases of nominal and zero LPOTP resistance. It 
is seen that the addition of the accumulators restores the system stability 
and is insensitive to reduced IPOTP resistance. 

1.4.7 Accumulator Failure 

A limited study was undertaken to examine the sensitivity of the 
system stability to isolated accumulator failures. The objective of the study 
was to obtain an indication of the criticality of such failures. With the 
available pitch-plane model it was possible to analyze the following two 
failure conditions: 

a. tckilure of both accumulators on the two lower engines 

b. Failure of the accumulate r on the upper engine 

Both of these cases were treated. The analyses were carried out with the 
L2 and the shifeed L54 modes; stability calculations were limited to nominal 
pump gains and downcomer structural mode amplitudes. In the case of the 
L2 mode it was found that loss of the accumulator on the upper engine had 
negligible effect on the system stability, and loss of the accumulators on the 
two lower enginfs resulted in a reasonably small loss in the stability level. 
These effects 'ire illustrated in Figure 19 where the system damping for the 
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Figure 20. Stability of L>54 Mode Shifted to 
Higher Frequency with Upper 
Engine Accumulator Ineffective: 
Rockwell Nominal Pump 
Compliances 



two failure conditions are presented together with the results for all the 
accumulators effective, 

Ir view of the results of the ground vibration test on Orbiter Vehicle 
101 (see discussion in Section 1.4,6), the analyses for the L54 mode were 
undertaken at the higher structural mo<le frequencies (25 — 32 Hz). To con- 
struct a critical stability case the Rockwell nominal pump compliances were 
employed to bring the first inter-engine mode of the propulsion s^'stem 
(without accumulators) into this frequency range. Since the upper engine has 
the predominant motion in the L54 structural mode it was not surprising to 
find that loss of the accumulators on the two lower engines did not produce 
any significant degradation of the system stability, whereas loss of the upper 
engine accumulator precipitated an instability (see Figure 20). The propulsion 
system mode involved in this instability had a strong inter-engine mode 
character. 

The above results should not construed to imply that the accu- 
mulators on the two lower engines are unnecessary. Clearly, the absence of 
these accumulators can be expected to yield a much greater stability loss in 
the L2 mode for conditions of larger downcomer motion and higher pump gains 
(see Figures 13 and 14), In addition, the role of these accumulators may be 
significant in maintaining stability for coupling between inter-engine propul- 
sion modes in the yaw plane and local engine structural modes involving one 
or both lower engines. Study of this latter possibility will be undertaken 
when a three-dimensional model of the MPS and vehicle structure is available. 

1. 5 SUMMARY 

The stability analyses undertaken for the liftoff event with the model 
5. 3A structural modal data set revealed the following for the system without 
accumulators : 

a. Degradation of the L2 structural mode (2,54 Hz) stability 
relative to the previous predictions made with the Shuttle 
model 4, 3 modal data (see Figure 8) 
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b* The presence of an instability in the L54 local upper-engine 

structural mode (17.8 Hz) due to coupling with the first inter- 
engine mode of the propulsion system (see Figure 10) 

c. Sensitivity of the LZ mode stability to variations in pump 
gains and in downcomer amplitude (see Figures 13 and 14) 

d* Aggravation of the instability in the L54 mode with decreased 
resistance of the LPOTP (see Figure 16) 

3 3 

With the addition of the baseline 0.017 m (0.6 ft ) gox accumu- 
latorSy it was found that satisfactory stability existed in all situations inves- 
tigated. Specifically 


a. The stability of the LZ mode was restored to a satisfactory 
level (see Figure 11). 

b. The instability of the L54 mode was eliminated and the damping 
in this mode restored to a satisfactory level (see Figure IZ). 


c. 


d. 


Th ~ '"'sitivity of the LZ mode stability to variations in 
r dn and in downcomer amplitude was greatly reduced 

( ure 14). 

Th j ility of the L54 mode was insensitive to decreased j 

rcsista ce of the LPOTP (see Figure 16). 


Based upon the location of a 30 -Hz local upper-engine structural 
mode in the orbiter ground vibration tests, sensitivity studies were under- 
taken with the L54 mode shifted to about 30 Hz. The results indicated mar- 
ginal stability conditions for the system without accumulators if the first 
inter-engine mode frequency was also in this frequency range. However, 
addition of the baseline accumulators again restored the stability to a satis- 
factory level even in the absence of LPOTP resistance. 


In support of accumulator charging studies by NASA and Rockwell, 
an analysis was cond cted of the variation of LZ mode stability with accumu- 
lator volume. It was found that an accumulator compliance equivalent to a 

3 3 

gox volume in excess of 0,007 m (0. Z5 ft ) was required at liftoff to avoid a 
50 percent loss in system damping for the assumed structural damping ratio 
of 0.01. 
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The effect of failure of the accumulator on the upper engine or both 
accumulators on the lower engines was examined for the L2 and L*54 modes. 

It was found that such failures had relatively little impact on the stability of 
the L2 mode for nominal pump gains and downcomer amplitude. Any accu- 
mulator failure^ however, is expected to lead to significant L2 mode stability 
loss for higher pump gains and downcomer motion (see Figure 19). Failure 
of the upper-engine accumulator precipitated an instability in the L54 mode 
(see Figure 20) with the Rockwell nominal pump compliance values. 

Finally, it should be noted that the present analyses did not involve 
the separation mode of the propulsion system (see Section 1. 3). Previous 
investigations (Refs. 4 and 7) have demonstrated the benefit of accumulator 
resistance in the elimination of an instability associated with this particular 
fluid mode. The matter of accumulator resistance will again be addressed in 
the next phase of our pogo studies. 
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2. APPUCATION OF QUARTER-SCALE LPOTP 
TRANSFER FUNCTIONS 


Pump dynamics play a very important role in pogo stability analysis. 
The pump model employed in previous studies (Refs. 1 through 4) does 
represent the existing state of the art (Ref. 5). However, the validity of this 
model is questionable because it was not based upon detailed pump experi- 
mental data. A beginning of the process for refining the pump model is 
described in this section using the pump dynamic transfer functions deter- 
mined experimentally at the California Institute of Technology (Caltech) for a 
quarter-scale model of the LPOTP impeller (Ref. 6). 

2. 1 CALTECH EXPERIMENTAL PROGRAM 

An experimental program for the dynamic response of cavitating 
turbomachines has been under way at Caltech, under contract from the NASA 
Marshall Space Flight Center (Ref. 6). In this study, elements of the linear- 
ized four-terminal transfer matrix of a quarter- scale LPOTP impeller have 
' ! been measured at a number of operating conditions using water as the working 

fluid. Tabulated results were made available to us by Caltech, including 
some results not published in Ref. 6. The study reported here employs four 
sets of transfer matrix results, each corresponding to a particular value of 
the pump cavitation number (specifically, = 0. 513, 0. 114, 0.046, 

0.024), all for a flow coefficient of 0.07 at 9000 rpm. 

A schematic drawing of the Caltech pump test facility is shown in 
Figure 21. The pump, driven by an electric motor, is located at the lower 
left. An upstream smoothing section, containing the upstream pressure 
transducer, precedes the upstream flow measurement by a laser doppler 
velocimeter. Likewise, a downstream smoothing section, containing the 
downstream pressure transducer, precedes the downstream flow measure- 
ment. Also, an accelerometer is employed to measure the pipe motion in the 
direction of flow for each flow measurement. This enables the absolute flow 
measurement by the laser to be converted to relative flow for the transfer 
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Schematic Drawing ot Dynamic Pump Test Facility (from Reference 6) 




matrix determination. Upstream and downstream, siren-type valves are 
employed to create periodic flow disturbances, with their individual amplitude 
and their relative phase adjustable to achieve various independent states of 
excitation. Other details of the setup are described in Ref. 6. 

For each steady operating condition of the pump (rotating speed, 
flow coefficient, cavitation number), and for each of seven discrete frequen- 
cies of excitation, data is collected at three or more independent states of 
excitation for the transfer matrix determinations at those frequencies. The 
data consist of pressure, absolute flow, and Une acceleration at upstream 
and downstream positions relative to the pump. By a least- square fit pro- 
cedure described in Appendix D of Ref. 6, the real and imaginary parts of the 
four elements of the following transfer matrix are determined: 
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whe re 

P = a dimensionless sinusoidal pressure 
m = a dimensionless sinusoidal relative mass flow 

u = upstream 

d = downstream 

w = the dimensionless frequency 

The dimensionless quantities are defined by 


P = >Vl '’"up 

m = Q/A;U,.p (5) 
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volumetric flow 
fluid density 

tip speed of the pump inducer 
inlet flow area 
physical frequency 
tip spacing of the inducer 


Studies were initiated for two possible methods of using the experi- 
mental results to describe the LPOTP for stability analysis. The first 
method consists of obtaining fitted polynomials (in the complex frequency, s) 
to the measured transfer functions. The second method consists of obtaining 
a lumped parameter, physically based mo<lel which reflects the experimental 
results. In both cases the data are employed in dimensionless form. Con- 
version of the dimensionless parameters to dimensional ones is defined in 
Table Z. 

2 . 2 CURVE FITTED TRANSFER FUNCTIONS 

This section describes the first attempt to employ the quarter-scale 
LPOTP results for stability analysis. The approach consisted of obtaining 
analytical fits of the measured transfer functions and using the resulting 
expressions to describe the LPOTP in the stability analysis equations. Two 
fitting approaches were tried: (P global polynomial fits valid for the entire 

frequency range and (2) piecewise linear fits valid over a limited frequency 
range. These direct analytical fitting approaches had the advantage of 
enabling a relatively rapid incorporation of the experimental data into the 
analysis. It was, however, somew \at crude and, when implemented, was 
judged to be of questionable validity. Nevertheless, the approach is outlined 
and the concerns are identified as background for possible future studies. 
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Table 2. Conversion Factors to Obtain Dimensional 
Quantities for Pump Modeling 


Dimensionless 

Quantity 

Multiply by Factor Below to Obtain 
Dimensional Quantity* 


SI Ur. ts 

Engineering Units 

Frequency, W 

U^^y2nH = 96. 1 3 Hz 

(same) 

Resistance, R 

PU^.p/2A. a 4.33 MN s/m^ 

Utip/^gAi * 0.285 sec/in.^ 

Inertance, L 

PH/2A. = 7. 16x10"^ MN s^m^ 

H/?.gA. a 4. 73x10"'^ sec^/in.^ 

Compliance, C 

2HA^PU^.p = 0.0924 m^/MN 

?.gHA./Uj.p = 5. HOxlO"^ in.^ 


Conversions based upon 

U.. = 3S. 9 m/s (1 17. 8 ft/sec) 

tip 

= 41.4 cm^ (6.41 in,^) 

H = 5. 94 cm (2. 34 in.) 

f* - 998 kg/m^ (Pg = 0.0361 Ib/in.^) for water at 21°C (70°F) 
= 9.81 m/s^ (386 in./sec^) 


g 




2 . 2 . 1 


Experimental Transfe** Functions and Global 
Polynomial Fits 


> 


( 


! 

*J 


The experimental transfer functions defined by Eq. (4) were adjusted 
by Caltech to remove the contrib\ition of the facility portions between the 
measuring stations and the pump itself. To distinguish these transfer func- 
tions from the nonadjusted ones, we use the subscript p on the Z's to denote 
pump: 




• 



l+ZpllM 


Zp2l'<sS 






1 

k J 


( 6 ) 


For our purposes the form was altered to obtain a correspondence with the 
form of the elementary model described by Eq. (1). Namely 


Pd = 


1 + z -P21-Pl2 


p - t 

■-pi 2 

u € 

.^■^-p22 . 


(7) 


Q, = {1+Z Q f £Z P 
d —pZZ u — p21 u 


where 


f = 



Note that the pressures and flows are now in dimensional form. 

The variations with frequency of the transfer function combinations 
appearing in Eq. (7) are shown in Figures 22 through 28 for the case of a 
cavitation number (J^ equal to 0. 114 (a repreaentative value for flight condi- 
tions), The figures show the experimental data points, a global polynomial 
fit of the data points, and, where applicable, the corresponding variation 
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from the elementary pump model used in past analyses. The global poly- 
nomial fits are for s = and have the form 


n 


Real part = c. 

j=0 


n 


Imaginary part = ^ ^ Cj s*^ 

j=l 


(j even) 


(j odd) 


( 8 ) 


'J 


The value of n employed is stated in each case to define the level of the least 
squares fit. 

Figures 22 and 23 show the real and imaginary parts of the transfer 
function ^^^2' dashed line in the figure represents the data fit to the order 
of n give.i on each figure. This transfer function has been referred to as a 
'*mas8 flow gain factor” (Ref. 13), and there is no equivalent term in the 
elementary pump model (i.e., Z .^^—►0). Figures 24 and 25 show the real 
and imaginary parts of the function dashed line in Figure 24 repre- 

sents an n = 2 fit of the real part data, the dashed line in Figure 25 represents 
an n = 1 fit of the imaginary part data. Figure 25 also contains the variation 
of the correspc ^mg term in the elementary model ^p2 
previously employed value of compliance for this cavitation number. The 
figure shows that the elementary model exhibits the same general trend as 
the experimental data. However, comparison shows that the value of cavi- 
tation compliance employed for the elementary model was considerably 
greater than that exhibited by the experimental data. This comparison could 
be influenced by several unknown factors so no conclusion is drawn. First, 
it has been shown theoretically that what has been called cavitation compliance 
for the elementary model includes a contribution of the mass flow gain fac- 
tor (Ref. 13). On this basis the qualitative relationship seen in Figure 25 is 
to be expected. Other unknowns are the effect of lox versus water for the 
working fluid and the effect of dimensional scaling. 
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Figure 24. Polynomial Fit (n = 2) of Real Z 



MODEL FREQUENCY. Hz 


Figure 25. Linear Fit (n = 1) of Ima , inary z 



1 } 


Figures 2b and 27 show th*» real and imaginary parts of pZ2^ 

an n = 6 fit of the real part data is shown in Figure 26; an n = 1 fit of the 
imaginary part data is shown in Figure 27. The variation of the corresponding 
terms employed for the elementary model is also shown. In the case of the 
real part, the elementary model term is simply the negative of the pump re- 
sistance; and in the case of the imaginary part, the term is the negative of the 
product of fret e^'cy and the pump inertance. On Figur? 26 it is seen that the 
constant oump resistance employed for the elementary model does not reflect 
the frequency dependence apparent in the experimental data. The most 
dramatic aspect of this dependence is the rapid decrease of the real part at 
the higher frequencies. As wa • noted earlier (Section 1.4.5), this behavior 
led to the performance of some stability analyses (for a high-frequency struc- 
tural mode) with an assumed LPOTP resistance of zero. The indication from 
Figure 27 is that the inertance associated with the experimental data is con- 
siderably larger than the value employed for the elementary model. 

Finally, Figure 28 illustrates both the real and imaginary parts of 
the function 1 + * -?p 21 — pl 2 /^^^— p 22 ^* equivalent elementary pump 

model term is the gain pump m + 1 , which should be compared to the real 
part data. For this particular test condition, the data indicates a pump gain of 
about 0. 9. Data obtained at lower values of the cavitation number have 
exhibited values as large as 2 , 8 . 

2.2.2 Piecewise Linear Fits 


In addition to the polynomial type of global fit, piecewise linear fits 
were obtained over a limited frequency range for s = iw. The form employed 
was 


Real part = + c^(s| 

(9) 

Imaginary part = + ^ 3 ®/^ 

where the coefficients c^ are real quantities. There is no point in showing 
these fits since they are essentially coincident will the data in the limited 
frequency range of applicability. 
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N.ODEL FREQUENCY, Hz 


Figure 26 


Polynomial Fit (n = 6) of Real 

tZpu'<**2p22)> 


Figure 27 


Linear Fit (n = 1) of Imaginary 

I^pl2''<‘"2p22>l 


-pll -pl2 -p21 


ELEMENTARY MODEL (PUMP GAIN) 
(REAL ONLY) 



0 20 40 


MODEL FREQUENCY. Hz 

Figure 28. Polynomial Fit of [ 1 +Zpj 2 - 
-ol2^^^‘'’-p22^^ 


IMAGINARY PART 



Figure 29. Propulsion System Eigenvalues Cal 
culated with Piecewise Linear Fits 
of Caltech Transfer Functions 


I 


Propulsion System Modes 


Concerns about the functional fitting approach arose when it was used 
to represent the LPOTP and then the modes of the propulsion system were 
calculated. When these modes were calculated with either the piecewise linear 
representation or with the global polynomial rep escntation, it was found that 
unstable modes were predicted for the propulsion system. Although instability 
is physically possible since the pumps are active elements, the results were 
viewed with suspicion. One prediction, utilizing a sequential piecewise linear 
fitting of the experimental data, is shown in Figure 29. The eigenvalues of 
the propulsion system are plotted in the complex s plane. The very strong 
instability (^-—0.7) was judged to be unreasonable. 

A concern with any fitting approach stems from the usage of data 
derived under simple harmonic conditions (s = ic«i) to apply to transient condi- 
tions (s = + iM. This was emphasized further by the experience gained with 

the polynomial fitting approach. In this instance it was found that the use of 
higher-order polynomials (say, n ^ 6 or 7) to fit the transfer functions re- 
sulted in the prediction of propulsion system instabilities at low frequencies 
even though the data fits for harmonic conditions looked completely reasonable. 
In retrospect, the difficulty may have been due to the fact that flow continuity 
through the pump was not constrained to exist as the frequency went to zero 
[that is, see Eq. (7)]. Thus, the fits should have been con- 


strained such that 1+Z 

— p22 


►1 and Zp2l^^ Figures 22 and 24 show 


that the real parts of the fitted polynomials did not satisfy the continuity 
requirement. This matter will be explored further in future studies. 

2.3 LUMPED PARAMETER MODEL 


2.3,1 Description of Model 

The lumped parameter model (in electrical analogy form) shown in 
Figure 30 was developed based upon an initial assessment of the test data. 
The model of the pump (LPOTP) differs from the usual elementary model of 
Eq* (1) only in that a parallel resistance-inertance circuit (resistances 2R, 
inertances L and aL) is employed for the pump flow impedance. The other 
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two pump elements are the usual cavitation compliance C and pump gain 

P 

m+1. The upstream facility section (containing a flow smoothing section with 
its pressure transducer P and an inlet feedline with its flowmeter, m ) is 
represented by the inertance and compliance C^. It is assumed for the 
model that the upstream state variables P^, m^ are physically located at the 
smoother /feed line flange interface (see Figure 21). The physical separation 
of the actual pressure and flow sensors is believed to be of second-order 
importance as was assumed in Ref. 6. The inclusion of the compliance 
reprebcnts to first order the small compliance of the smoother transition 
section and the feedline; this compliance is to be determined from the test 
data, fh: upstream inertance = 10. 7 is taken to be that given by the geo- 
metry between the smoother/ feedline flange and the beginning of the pump 
blades. The downstream facility section is modeled in a corresponding fashion 
with inertance = 4. 1 based on geometry and the compliance to be deter- 
mined from the test data. The compliance is relatively large since it 
includes the effective compressibility of the relatively voluminous downstream 
smoothing section. 

The elements for the model of Figure 30, as defined in Eq. (8), 
are as follows: 

1 + Z,, = (m+l)B - w^(C +C B)(L + L .) + iWR (C +C B) (10) 

11 p u p d' p' p u ' 

Zj 2 = - RpB - iw[(m+l)L^ + (Lp+L^)B] (11) 

Z^i = «2(Cp+C^B)CjRp + iCp+C^B)(l . L^)l| 

(12) 

1 + Z 22 = [1 - ^^^C^(Lp+ L^)]B -w^(m+l)CjL^ + iwC^RpB (13) 
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Figure 30. First-cut Pump Mathematical Model 


NORMAIIZED REAl AND IMACINAkY PARTS 



n-wL<R 


Figure 31. Normalized Pump Impedance with 
Test Data under Noncavitating 
Conditions ( <Tp = 0.513) 





where 


^ ^ ^[2+g^d-fa^)] 


(14) 


uL 

E 

R 


ifl(l+a)(lfg^a) 
4 + fi^(l+a)^ 


(15) 


Q = cjL/ r 


B = 1 - U C L 
P u 


As previously noted, the values of and are found from feedline geo- 
metry* The values of the other parameters (C^, C, R, L, a, m+1) are 

to be determined from a matching of the above Z.^ elements with their experi- 
mental values over the range of test excitation frequencies. It is assumed 
that only the pump cavitation compliance and the pump gain m+1 vary with 
the pump cavitation nunme." As previously stated, data for four values of 

cavitation number are employed for this sfcv v. 

The basis for the parallel representation of the pump flow impedance 

Z *-> experimentally observed nature of the variation with frequency of 

pi Z 

its real and imaginary parts. From Eq. (11) we see that the impedance Z 


pl2 


is given for the fully wetted or noncavitating (C^ = 0, m + 1 = 1) conditions by 


pl2 


12 


i(*/( Li + Li , ) 

u d' 


(16) 


Figure 31 shows normalized Z dala for the fully wetted case = 0.513) 

pic P 

superimposed on a family of curves for the parallel circuit shown with a as 
the parameter. Two normalization factors are employed. One is the steady- 
state pump resistance R obtained in the usual fashion from pump steady per- 
formance data as the negative of the partial derivative of the pump head with 
respect to pump flow; from Figure 5. 2 of Ref. 6 the value R = 16. 6 is obtained. 
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The second normalization factor is r; = 0. 8 to account tor a suspected error 

in the calibration of the upstream flow measurement rn^^. This error was 

rationalized on the basis that Ihio flow was the most suspect measurement and 

that it is necessary for Re as so that inflow equals outflow under 

steady flow conditions. The experimental data indicated that Re — 0.2 

instead. Let us assume that the actual flow m = nm’ , where m' is the 

measured flow; the value 77 = 0,8 provides the indicated correction. Since 

the Z.. are determined using the measured flow rh‘ , the corrected transfer 
ij ^ u" 

matrix becomes 


' - 

Pd 

► 

'uz„ 


i 

f ^ 

p 

u 

,”d. 


/21 

d+Zjz'/". 


m 


(17) 


As a consequence, the experimental data for ^vere divided by * 0. 8 
before evaluated using Eq. (16). 

The value of a (between 0 and 1) determines the high frequency 
asymptote of the real part at a value of between one and two times the steady- 
state resistance R; this matches the trend of the test data shown on Figure 31 
and also the trend of data obtained in another test program (for example, 
see Figure 9 of Ref. 14). It is informative to note that a = 1 yields the expected 
result Zpj 2 - R + i(jL/2. The match between the data and the mathematical 
model shown on Figure 31 was achieved by trial and error using the value 
L = 320 which establishes the scaling of abscissa values. The correlation, 
most particularly for the imaginary part data, suggests the value a = 0.04, 

Note that these are the values appearing on Figure 30, 

Thus far, the values of L^ and were determined by feedline geo- 
metry, K from the pump performance data, L and a from the matching to the 
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circuit shown in Figure 31. The remaining parameters were evaluated as 
follows: 


m +1 = Re(l+Zjj) (low w) 


(18) 


ImZ__ 

C = 

— d (,/Re Z 


(low u) 


-12 


C 

— u 


+ C 
-P 


ReZ^l 


wlmZ22 


(JReZ 


(low u) 


(19) 


(20) 


For the fully wetted case (C^ = 0, m + 1 = 1), Eq. (20) provides the basis for 
the value of C^. Thus, all of the model parameters are determined. The 
values appear on Figure 30. 

2. 3, 2 Evaluation of Results 

Figures 32 through 39 show the calculated real and imaginary parts 
of the elements versus frequency frc^m the mathematical model in 
Figure 30. The results are encon . in that the model does qualitatively 

show most of the major trends in tne data, particularly for the two higher 
cavitation numbers. 

Available information on steady engine operation in flight yields a 
range of cavitation numbers from about 0. 1 to 0.4. Thus, the test condition 
cavitation number <7^ = 0.114 is the most representative of flight operating 
conditirns. A review of Figures 32 through 3 9 shows that the match for this 
value rf ik? relatively good. The matches for ImZjj and ReZ 2 j are par- 
ticula**'/ The matches for ReZ^ 2 > ^^—22" ^^—22 

bettei than tr:oy appear at first glance when the factor 1. 25 (equal to I/ 77 ) is 
applied to the data [Eq. (17)]. The other two cases, ReZj^ and ImZ^j, show 
moderately good correlation; in each case a constant factor would move the 
experimental data points into a much improved match to the mathematical 
model. 


2-lh 




2-19 






DIMENSIONLESS FREOUENCY.y 


Figure 34. Real Part of Z 
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DIMENSIONLESS FREQUENCY, w 


Figure 35. Imaginary Part of Z 
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For the most part the model provides a good match at low fre- 
quencies. Special attention is called to the initial shaj of ImZ^j which 
reflects the total compliance and pump gain as follows [^ee Eq, (12)]: 


^ ► (ni + 1) 


CJ- 


+ C 


c 

— u 


( 21 ) 


On the other hand, the poor quantitative high-frequency correlation 
at the lower cavitation numbers strongly suggests that further model develop- 
ment is needed for an adequate description of pump dynamics under conditions 
of major cavitation. For example, such conditions may occur during the 
after-SRB-separation transient and <Iuring a lox depletion where the lox sur- 
face moves well down the feedline prior to engine shutdown. 


2.4 SUM^ ^ \RY 

Two preliminary studies were made for the use of the Caltech ex- 
perimental results for mr- ’eling the LPOTP, Each study showed questionable 
results. Further efforts along both lines are planned for the future. 


In the first study two forms of functional fitting, global polynoraial 
fits and piecewise linear fits, were applied to the Caltech pump data. Con- 
cerns about the validity of the results for use in stability analyses arose 
because various instabilities were predicted for the main propulsion .system 
when the fitted functions were employed for the LPOTP. 


The second study was aimed at developing a lumped parameter, 
physicallv based model. The result was the model shown on Figure 30, A 
comparison of the results from this model with the measured data appears 
on Figures 31 through 3 9. The comjjarison shows an encouraging similarity 
in the major trends, particularly for the cavitation number most represen- 
tative of steady flight operation {(T^ - 0, 114) and higher. However, several 
matters need to be investigated and resolved before such a model can be 
accepted; 


a. The need for the correction of the upstream flow data to 
achieve continuity for scoady flow 
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b. The poor quantitative correlation at high frequencies at 
low cavitation numbers 

c. The use of the model to calculate hydraulic modes of the 
Caltech test system to verify a lack of major discrepancy 
with observed self-excited system oscillations. 
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3. INITIAL POC.O STUDY OF OMS OPERATION 


An exploratory evaluation of pogo stability during operation of the 
orbital maneuvering subsystem (OMS) is described in this section. The 
concern is prompted by the fact that two other vehicles powered by generally 
similar pressure-fed engines (see Table 3), namely the Delta Stage II and 
the French Diamant B Stage I, hav'e experienced strong pogo instability in 
flight (Refs, 15 and 16). The Delta Stage II, powered by a derivative lunar 
module descent engine (LMDE), has experienced sustained peak axial vibra- 
tions of up to 45 g at 125 Hz. The vehicle structural mode was a local one 
involving primarily engine axial motion with some interaction with the vehicle 
tanks. The Diamant B had experienced vibrations as high as 20 g at 45 Hz in 
the second longitudinal mode of the overall vehicle. The Delta Stage II pogo 
has apparently been corrected by ballasting the engine with a 19. 5 kr; mass 
(43 lb) to reduce the structural gain. The Diamant B instability was sup- 
pressed by a combination of stiffening of the junction of the thrust structure 
with the aft propellant tank (to reduce structural gain) and by bonding visco- 
elastic material around the aft half of the same tank (to increase structural 
damp'ng). It is noteworthy also that both the Delta Stage II and the Diamant 
B Stage I underwent pogo instability during stage static firings. The overall 
lesson for the Shuttle is that pogo prevention should be vigorously pursued on 
the OMS. 

The following sections deal with preliminary structural and propul- 
sion dynamic models, studie • to investigate mechanisms for potential insta- 
bilities, and an initial look at suppression for an identified potentiallv strong 
destabilizing mechenism. Details of the models are found in Appendices C 
and D. 

3. 1 QMS STRUCTURAL DYNAMICS MODEL 

The OMS pod illustrated in Figure 40 is assumed to be fixed at the 
orbiter fuselage interface with the major portion of its mass concentrated in 
the lMS and RCS propellant tanks, the OMS high-pressure helium storage 
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Table 3, Selective Comparison of OMS with Pressure-Fed 
Vehicles Experiencing Pogo F^enomena 



OMS 

Delta Stage 11 

Diamant B Stage 1 

Engine Mass 

113-122 kg 
(250-270 Ibm) 

104-132 kg 
(230-290 Ibm) 

(Unknown) 

Thrust 

27 kN 

44 kN 

40 kN 


(6000 Ibm) 

(10,000 Ibm) 

(9000 Ibm) 

Props Hants 

MMH, 

UDMH, 

UDMH. 


3070 Ns/kr, 

2960 Ns/kg 

2880 Ns/kg 

sp 

(313 sec} 

(302 sec) 

(294 sec) 

Frequency 

45 Hz 

125 Hz 

45-47 Hz 

Structural Mode 

Local Engine 

Engine with Tanks 

Second Longitudinal^ 
Large Engine Motion 








bottle, the RCS thruster arrangement and the OMS engine. The pod structure 
consists of aluminum and titanium truss members of skin and stringer con- 
struction, and reinforced panels. Mass and flexibility matrix data were 
based on a 123-DOF mathematical model supplied by the Rockwell International 
Space Division (RI/SD). 

A set of OMS pod structural dynamic rrodes (corresponding to the 
100 percent fill condition in all propellant tanks) was calculated by RI/SD. 

The tanks were assumed rigid with flexible aft bulkhead supporting structure. 
The lowest tank-dominated OMS longitudinal modal frequencies were at 10. 0 
and 12.6 Hz. In addition, a high-gain localized engine axial mode was cal- 
culated at 43. 5 Hz. 

3,1.1 Simplified OMS Pod Structural Modes 

An independent calculation of pod structural modes was carried out 
since the RI/SD supplied data were limited to the 100 percent fill condition 
and the assignment of a single pod structural mass point at the pod center of 
gravity was not thought to be adequate for pogo analysis. 

The RI/SD assumption of rigid propellant tanks was first examined 
with an elementary flexible tank model (Ref. 17) supported by a ring at the 
aft dome/cylinder interface which connects to a support spring. The modes 
of this simplified model verified that a rigid tank model is valid for all pro- 
pellant fill conditions for frequencies below 100 Hz, 

For the present investigation, a simplified 11-DOF longitu linal 
dynamic model with engine lateral dynamics wa sembled. The model with 
selected significant mass points is illustrated in , ,gure 41 and the RI/SD 
based mass and flexibility matrices are presented in Appendix C, In the 
absence of distributed structural mass data, a conservative approach was 
taken whereby the structural mass (lumped at its center of gravity by RI/SD) 
was completely ignored. The burn time dependent tank masses were taken as 
functions of propellant level with and V denoting the OMS and RCS fractional 
fill levels , respectively. 
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MODEL INCLUDES VARYING PROPELLAMT MASS 

OMS Fua S 2041 kg (4503 Ibm) RCS FUa 1 371 kg (817 Ibm) 

OMS OX < 3372 kg (7433 Ibm) RCS OX < 616 kg (1357 Ibm) 


Figure 41. 11-DOF OMS Structural Dynamic Model 



Figure 42, OMS Pod Structural Mode Frequencies vs. 

OMS Propellant Level (RCS 100% Full) 
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Calculation of the modes as a function of OMS propellant level was 
carried out. The mode spectra are presented in Figures 42 through 44 for 
three RCS propellant fill levels of 100, 50, and 1 percent, respectively. In 
general, the results indicate a persistancc of the localized high-gain engine 
axial mode at 45 Hz (versus RI/SD calculated 43. 5 Hz) for nearly all propel- 
lant levels. Some splitting with OMS tank-dominated modes occurs near tank 
depletion. Typically the modal mass associated with the engine localized 
mode is about 128 kg (282 Ibm) with ^ s 1 (slightly above the assumed engine 
mass of 113 kg (250 Ibm), Tank -dominated modes, which increase in fre- 
quency with propellant consumption, potentially contribute to pogo destabili- 
zation trends. In these modes engine gimbal motion is of the same order as 
the tank axial motion and high tank bottom modal pressures and motions pro- 
vide for strong propulsion system coupling. A summary of engine dominated 
(45 Hz) and significant tank dominated mode spectrum parameters are pre- 
sented in Appendix C. 


Conservative Upper Bounds on Modal Gain 


In view of the preliminary nature of the current OMS structural 
dynamic model, a simplified description resulting in a conservative estimate 
of modal gain was constructed for parametric study of the pogo phenomenon. 
If one assumes that the contributions to modal mass (i.e,, modal kinetic 
energy) are due only to engine axial motion and OMS propellant tank motion, 
the modal mass is expressed as 


(M„_ 


where the engine modal displacement is taken as unity {<f>^ * 1). The above 
expression represents a lower bound on modal mass for a given mode shape 
since typically there would be contributions due to other components (i.e.. 


RCS tankage). The lower bound modal mass yields an upper bound in engine 
modal gain (0^/M^U With the tanks contributing to the increase in modal 


mass (or decrease in engine modal gain), there are associated modal tank 
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QMS PROPELLANT FILL LEVEL 


Figure 43. OMS Pod Structural Mode Frequencies vs. 
OMS Propellant Level (RCS 50% Full) 
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OMS PROPELLANT FILL liVEL ^ 


Figure 44. OMS Pod Structural Mode Frequencies vs, 
OMS Propellant Level (RCS 1% Full) 
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bottom pressures (for closed tanks) which are expressed approximately as 

= wMAV P.o ' 

for the fuel and oxidizer tanks, respectively. Note that the tank length, 
is the same for fuel and oxidizer tanks. 

For the purposes of a preliminary sensitivity study, the following 
modal parameters are varied: (1) tank modal displacements and (2) modal 

frequency in the range of 20 to 65 Hz (within which tank-dominated and engine* 
dominated modes occur for th<5 11 -DOF model described in the previous 
section). 

3.1.3 Structural Mode Representation for Pogo Analysis 

"Tie representation of modes for closed-loop pogo analysis employs 
adjustments for open tank/feedline fluid interfaces and the associated reac- 
tions (Ref 1). For convenience, volumetric outflow displacements rela- 
tive to dotted tank bottoms have been chosen as generalized flow coordinwt .s. 
The resulting dynamic equations associated with a single structural mode 
expressed in a convenient form for organization of the stability equations are 






\ 

A 



R 


tf 



F ' 

e 

V 


tf 


tf 


(24) 
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where F*=F +F.+ F • These equations are subjected to the modal 
e e ef eo 

constraint 



In the present analysis the outflow inertances and are neglected since 
they are generally irs gnificant (Ref. 1), and the outflow resistances 
and R^^ are included in the feedsystem model and thus not included again in 
the structural model. 

It should be noted that the tank bottom reaction terms, F^^ and 
the generali'^ed interaction forces, F^^ and F and the tank bottom pres- 
sures, and ultimately do not appear explicitly in the final closed- 

loop stability equations. Equal and opposite interface loads are applied to 
the feedsystem, and oupling of the pod and feedsystem substructures results 
in the automatic cancellation of the explicit interface loads; they appear 
implicitly as mass, damping and stiffness coupling terms in general. This 
is the same mathematical phenomenon that occurs in modal coupling and 
substructure analysir; it is guaranteed as a consequence of the Lagrange 
equations. 

3.2 PROPULSION SYSTEM MODEL 

The OMS includes interpod and PBK-pod feedlines as well as intra- 
pod interconnections to the RCS. The present analysis is limited to the basic 
propellant delivery configuration. For the purpose of the current early pogo 
study, the geometrically complex feedlines and engine manifolding are 
assumed subject to axial-only structural motion excitations in a manner 
consistent with the pod structural model. Moreover, since feedline bracketry 
was not finalized when configuration data were obtained, the excitations are 
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limited to tank support structure axial motion and engine gimbal axial motion. 
A lumped parameter, finite element description of the feedlines and engine 
manifolding was chosen due to the geometric comp’^xity of the system. 

The remaining component of the propulsion system model is the 
description of combustion dynamics. The linear combustion model includes 
constants derived from steady performance data and the usual time lag and 
delay constants. 

3.2.1 Feedsystem 

The intrapod feedlines for fuel and oxidizer consist of all line sec- 
tions between the tank interfaces and engine manifold interfaces including the 
side branches leading to closed crossfeed valves and line sections on the 
OMS engine. The three-dimensional geometry of the feedlines is approxi- 
mated by a planar, longitudinal- lateral configuration with the longitudinal 
runs properly scaled. Such an approximation provides for proper accounting 
of the longitudinal hydraulic acceleration heads required for the present 
idealization. 

The lumped parameter, finite element feedline and engine manifold 
models are illustrated in Figures 45 and 46 for the fuel and oxidizer circuits, 
respectively. Hydraulic data for the fuel and oxidizer circuit models and 
details of the models are presented in Appendix D. The chosen set of dyna- 
mic variables consists of (1) volumetric flow displacements Q at fluid nodes 
taken relative to local line structural motion and (2) line structural displace- 
ment X . Assuming generally negligible longitudinal flexibility of the line 
s 

structure and allowing rotation of cornc'rs due to restrained bellows (univer- 
sal joints), one finds tne feedline dynamic equations take the form 
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Figure 45. OMS Fuel Circuit Model Schematic 



Figure 46, OMS Oxidizer Circuit Model Schematic 
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The feedsystem models for fut;! and oxidizer circuits are assembled by a 
multistep process. Utilizing the fuel circuit (Figure 45) as an example, the 
assembly steps consist of: 

a. Formation of simple node chains, namely 1-6, 7-l6, 17-21 and 
22-25 with 26-28 reserved for the line structural motions 

b. Application of continuity constraints and support constraints, 
namely 

e,7 = e(,-e, 

' ®16 



c. Incorporation of Guyan reduction (Ref. 18) retaining DOF 1,4,6, 
7, 9, 11, 16, 23-25, 26-28, as indicated in Figure 45. 

All constraint and reduction transformations are applied in a symmetric 
manner as for typical structural dynamic models. 

In the present application, the generalized mass matrix partition 
corresponding to the dynamics of the line vith frozen fluid mass is 
neglected implying that it is accounted for in the pod structural model. The 
load distribution vectors and consist of single unit entries at rows cor- 
responding to the tank and chamber interfaces, respectively. The load 
vector Fg is comprised of the appropriate structural interface reactions, 
such as F^j and F^^ for the fuel circuit. 

3.2.2 Combustion Dynamics 

The well-known linear characterization of the relationship between 
chamber pressure and injector face propellant flow perturbations is of the 
form 

-Tj.s ^ s ^ 

(1 + T s)P = R e e.f+R e G. (28) 

reef if CO 10 

with s representing the Laplace variable. The chamber resistance coeffi- 
cients, and *re estimated on the basis of steady engine performance 
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data (Ref. Z), which in the present case was obtained from actual engine 
development tests. The chamber residence time is calculated from 
chamber geometry data and the combustion time delays, and have 

been provided by Aerojet Liquid Rocket Company. Our studies showed that 
the amplification and phase lead due to the Klystron or clumping effect 
(Ref. 19) was negligible and it was thus ignored. The values of all combus- 
tion constants along with some basic performance data are presented in Table 4. 

Since the present formulation, and the analysis approach discussed 
in the next section, require a set of first-order, real, constant-coefficient 
differential equations, an approximation of Eq. (28^ has been derived The 
approximation, first used by us during Delta Stage II pogo studies in 1975, 
employs fourth-order Fade expansions (Ref. 20) for the exponentials and a 
state variable transformation (Ref. 21). The resulting state equations are 
described in detail in Appendix D. 


3.?. 3 Propulsion System Modes and Frequency Response 

The modal and frequency response characteristics of the propulsion 
system were studied prior to actual closed -loop stability calculations to 
obtain a basic understanding of structure /propulsion system interaction 
mechanisms and data for preliminary open-loop stability analysis. To this 
end, the assembled fuel and oxidizer circuit models with excitation at the 
structural interfaces, Eq. (26), were coupled with the combustion dynamics 
representation, Eq. (28), 


(L)|©| 
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(29a) 


(29b) 


The matrix consists of of Eq. (26) augmented by zero columns cor- 
responding to the engine state variable s X which are X X of Eq. (D- 1 3). 

c 1 o 
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Table 4. Engine Dynamic Combustion Parameters 


Steady State 

Thrust = 27 kN (6000 Ib^) 

P =0.86 MPa (125 psia) 
c 

Wj, = 3.28 kg/sec (7.22 Ib/sec) 
clf 

W, = 5.40 kg/sec (11.91 Ib/sec) 
do 


=310 cm^ (48 in,^*) 

MR = 1.65 

I = 3070 N-s/kg (313 sec) 

sp 

— = 1.01 sec/cm^ (6.53 sec/in. 

A^g 


Perturbations 


d I 

s 


d MR 


% = 142 N-s/kg (14. 5 sec) 


]— 


❖ 

C R 7 

^ a 0.8 sec/cm'^ (5.21 sec/in.") 

A^g Pfg 

, — — = — — = 1, 14 sec/cm^ (7. 3 sec/in. ) 

L "o8 


Time Constants 

Residence Lag = 0*88 ms 
Fuel Delay « 0* 6 ms 
Ox Delay s 1. 0 ms 





! 


f 
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I 




r 


Transformation to a convenient s«'t of first-order dynamic equations is 
effected by definition of the flow velocity variables 

[v| = |e| ( 30 ) 

yielding after some manipulation 



The complex prooutsion system modes exhibit coupled dynamics of 
the fuel ane o-.idizer cir/ hae to bipropellant combustion dynamics. For 
basic undei standing of the ^^ynamics of the individual propellant circuits it is 
necessary to first examine the system modes for fuel ^nd oxidizer circuits, 
separately, in the absence of combustion (that is, = 0 with the lowest 
row block partition : £q. (31) removed). The fundamental complex fuel 

and oxidizer modes (pr edominantly real) are illustrated in Figures 47 and 4P 
In both cases the modes indicate a major contribution of the side branches 
to the dynamic behavior. T* 3 low^er propulsion system modes and eigenvalues 
with combustion dynamics included are summarized in Table 5 along with the 
corresponding modes without combustion. Althoxigh the fuel and oxidizer 
circuits are somewhat coupled in these modes, the complex eigenvalues and 
mode shapes are not greatly affected by combustion dynamics. This insensi- 
tivity v/as due to the relativ'ely low impedance associated with the combustion 
parameters; in parti<*ular, the chamber resistances are significantly lower 
than the propellant circuit ilow resistances localized near the injector. The 
ratios are 




R /R S 0.56 
CO ox 


' 32 ) 
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Table 5. Summary of Lower PropuLion System Modes 


Without Combustion (P = 0) 
c 


Eigenva lue 


With Combustion ( using Eq. (28)] 


Eigenvalue 


a je. 

if lO 


Fuel rigid 


Ox rigid 


8^0 


Ox overdamped 
Fuel overdamped 


Ox fundamental 


Fuel fundamental 


Ox second 


Fuel second 


8 * -94,02 
8 « -176. 5 


Overdamped 

Overdamped 


( = 0,101, f a 67,6 Hz Ox fundament il 0,27, -130 deg 

i 9 0.113, f 79,6 Hz Fuel fundam*enta 1 16,7, 1 9 deg 

( = 0,105, f = 125, 9 Hz Ox second 0, 58, -132 deg 

( - 0,092, f a 134,0 Hz Fuel second 3.30, -25 deg 





The frequency response of the propulsion system subjected to 

individual unit structural excitations, X , was calculated utilizing t!ie com- 

s 

plete set of p. opulsion system modes (by the method described in Appendix E), 
The chamber pressure response to fuel tank, oxidizer tank, and engine axial 
excitation are illustrated in Figui*es 49 through 51. Ail responses 
exhibit fairly flat gain variation below o5 Hz with phase in the range of -30 
to -140 deg. The responses to tank acceleration are shown for tank propel- 
lant levels of zero and 76. 2 cm (30 in.). The effect of tank bottom pressure 
is accounted for in the frequency response by augmenting the acceleration 
head terms with the appropriate tank liquid heights 

htf = (33) 

noting Eq, (23). On the basis of the nearly constant difference between the 
responses with and without tank fluid, it is estimated that the effective upper 
line acceleration heads are approximately 71. 1 cm (28 in.) for both the fuel 
and oxidizer lines (h^^ and h^^, respectively). 

3.3 STABILITY ANALYSIS 

The present stability analysis includes: 

a. An initial open- loop investigation for identification of potential 
destabilization mechanisms 

b. A closed-loop analysis illustrating the nominal stability 
associated with the engine localized 45-Hz structural mode 

c. Worst-c^ "5e damping loss envelopes 

All of the above employ nominal propulsion system parameters. 

3.3.1 Open- Loop Stability Approximation 

Approximate open-loop stability analysis provides basic under- 
standing of the pogo feedback ^chanism and the role of system parameters. 

In the present study, open-loop analysis is extensively used in nominal para- 
meter and sensitivity calculations and in a preliminary worst-case investi- 
gation. 
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Figure 49. Chamber Pressure Response to Fie, «re 50. Chamber Pressure Response to Oxidizer 

Fuel Tank Motion (dP^/aX^^) Tank Motion (dP^/dX^^) 



Figure 51. Chamber Pressure Response to Engine 
Gimbal Motion (^P /d^ ) 



I 


The approximate equation for response of a single structural mode, 
if one ignores all interaction forces other than thrust perturbation, is 


X t CJ X + 
e o o e o e 



(34) 


with engine axial motion as the generalized displacement (mode normalized at 
X^) and representing the thrust coefficient given in Table 4. The total 
chamber pressure response due to structural mode excitation is 


P 

c 


bP 



e 





'Vo / 


X 

e 


(35) 


for which the component response functions are presented in Figures 49 
through 51 (for h^^ = h^^ = 0), If Eq. (3 5) is substituted into Eq. (34) at the 
structu»*al resonance U ^ a simplified estimate of closed-loop damping is 
obtained 




(36) 


based on the open-loop damping loss 








+ 



(37) 


with the asterisk ("i^) denoting the response component at -90 deg phase. As 
mentioned earlier, it should be noted that the above open-loop stability 
approximation accounts for the thrust feedback interaction loads applied to 
the structure alone. Reaction forces associated with open tank bottoms and 
feedsystem pressures at bends are ignored. A consistent account of all 
interactions inherent in the component models is realized in the closed-loop 
procedure outlined oelow and discussed in Appendix D. 
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3.3.2 Closod-Loop Formulation 

The closed-loop stability equations are constructed in a manner 
similar to the response equations of Section 3.2.3. The procedure consists 
of (1) assembly of the combined OMS pod/feedsystem hydroelastic equations 
and (2) assembly of the first-order coupled hydroelastic/combustion dyna- 
mic equations. 

Assembly of the OMS po<l/ feed system hydroelastic equations is 
accomplished by enforcement of modal and boundary constraints on the OMS 
pod structure, Eq. (24), and on the feedsystem, Eq. (26). Details of the 
constraint matrix as formed in the computer calculations are presented in 
Appendix D, Upon application of the constraints to Eqs. (24) and (26) in a 
symmetric manner, the following consistent set of hydroelastic equations is 
obtained; 
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(38) 
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The mass coupling partition (/?h)' (A represents the generalized acceleration 

head contributions due to feedsystem and tank heads combined, and the 

generalized thrust chamber loading matrix P' consists of non-zero entries 

se 

Pgg (chamber pressure to fuel circuit) 

P' (20, 1) a -1. (chamoer pressure to oxidizer circuit) (39) 

s c 

P' (21, 1) = A (thrust) 
se e 

Due to the consistency of constraints, all other explicit interface 
reaction terms do not appear in Eq. (38) (i.e., equal and opposite interactions 
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The closed-loop stability equations consisting of Eqs, (38) and (29b) 
(with slightly modified du<‘ to one additional variable) are finally constructed 
by introduction of the generali^^ed velocity vector 



(40) 


and some matrix manipulation resulting in 



(41) 


with and representing the generalized mass, damping and stiffness 

matrices defined in Eq, (38), The closed-loop matrix equation set, Eq. (41)^ 
is in the form of a standard algebraic eigenvalue problem and the solution 
yields the complete set of complex eigenvalues and eigenvectors. 


3,3.3 Nominal Stability Analysis Results 

Open-looo stability estimates were first calculated for the primary 
OMS pod axial mode (with RCS 1 percent full) and the engine axial gimbal mode 
(with RCS 50 percent full) versus OMS lank fill condition. The basic data 
for these calculations consists of modal data presented in Table C-3 (along 
with a fraction of critical viscous damping = 0.01) and the propulsion 
system response data presented in Figures 49 through 51. The resulting 
nominal stability versus OMS tank fill level, illustrated in Figure 52, indi- 
cates a nearly inva’^iant ermine gimbal mode instability due to the localized 
nature of the 45 mode and a destabilizing trend in the primary axial mode 
with a minimum 8. 9 dB margin near the 10 percent OMS fill lev^el. The dB 
margin is called the damping gain margin and defined in Eq. ( 3 ) of Ref. 5. 

margin = 20 log z ( 4” < 4” > ^42) 

lOA-A ’c’o 

‘»o ’c 
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Figure 52. Nominal Stability vs. OMS Tank Fill Level 



Figure 53. Stability Sensitivity of Engine Gimbal Mode to 
Frequency Variation 
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In the engine gimbal mode, th<^ thrust feedback primarily due to dP^/dX^, 
has a nominal phase of 214 deg lagging the stabilizing structural mode 
damping force, Tht* stability of the primary axial mode is de- 

graded by major components of all three chamber pressure response func- 
tions which produce a net feedback opposing the structural mode damping 
force. An instability is not realized due to the low structural gain associated 
with this mode. 

limited closed-locp stability calculations were performed for the 
primary axial and engine gin.bal modes, each at particular tank fill levels. 

The closed-loop damping associated with the engine gimbal mode (Figure 52) 
is extremely close to the open-loop estimate since interaction forces 
neglected in the open- loop approximation are insignificant. On the other 
hand, in the primary axial mode, the open- and closed-loop damping esti- 
mates differ somewhat due to neglected interaction forces. 

A detailed look at the closed- loop eigenvector (and similarly the 
propulsion system frequency response) indicates that the 45-Hz instability 
and the degraded stability of the primary axial mode are bipropellant phe- 
nomena, In particular, for the 45-Hz mode, the fuel to oxidizer injector 
flow ratio about 0. 9. For the primary axial mode (10 percent 

OMS fill, 1 percent RCS fill), the flow ratio is about 0.8, This observation 
suggest.^ the necessity of accumulators for both fuel and oxidizer circuits 
for adequate suppression, 

3,3.4 Variation of Engine Gimbal Mode Frequency 

The sensitivity of the engine gimbal mode to medal frequency varia- 
tion was examined on the basis of the open-loop estimate. Modal parameters 
associated with the OMS tank 75 percent fill level and RCS 50 percent fill 
level were used. The results presented in Figure 53 illustrate that instability 
may occur over a wide range of gimbal frequency (15-50 Hz), Moreover, 
for the present nominal set of propulsion system parameters, the system 
becomes phase stable as engine gimbal frequency approaches the fundamental 
oxidizer mode at 67. 6 Hz. It should be emphasized that this higher frequency 
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trend does not necessarily suggest stiffening of the gimbal structure for pogo 
alleviation. Variation of combustion time delays and/or propellant bulk 
moduli off the current preliminary values can effect unfavorable phase shifts. 

3.3.5 Damping Loss Envelopes 

In view of the preliminary nature of the OMS pogo model, conserva- 
tive worst-case damping loss envelopes are useful for definition of suppressor 
requirements. On the basis of the conservative upper bounds on modal gain, 

Eq. (Z2), and the representative upper-bound (phase varied) propulsion sys- 
tem response functions, an open- loop based damping loss envelope is defined as 




dP 


dX 


max 


*«i (‘ ^ 








max 




dp 


dx 


to 


max 


M 




(43) 


The chosen upper-bound propulsion system response values, valid in the 
15-60 Hz frequency range, are 


dP 

t 

dX 


= 110 Pal {ml s^) or 4x10”^ psi/(in. Isec^) 


max 


dP 


dX 


tf 


54 Pa/(m/s^) or 2x10 ^ psi/ (in./ sec^) 


max 


(44) 


dP 


dX 


to 


= 220 Pa/(m/s^) or 8x10 psi/(in,/sec^) 


max 


^-24 


Further simplification of the damping loss envelope expression is accom- 
plished by separation into somewhat more conservative individual components 
due to engine motion, fuel tank motion and oxidizer tank motion, respectively, 

i. e. 



In the case of damping loss due to engine motion, a reasonable lower-bound 
modal mass was chosen to include fuel tank inertia only; it should be noted 
that the fuel tank is always lighter than the oxidizer tank. 

Envelopes of the individual damping loss components are presented 
in Figures 54 through 56 as functions of modal tank motion with OMS tank fill 
condition as a parameter. The damping loss due to engine motion, 
decreases with increasing tank motion and fill level as expected. The damp- 
ing ratio losses due to tank motion, and increase proportional to 

<t>^ (tank modal amplitude) for low values of 0^ and asymptotically decrease 
inversely proportional to 0^ for large values of 0^. The magnitude of 
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Figure 54. Envelope of Damping Loss due to Engine Gimbal 
Motion vs. Tank Modal Amplitude 








greater than primarily due to the higher density of the oxidizer which 

manifests itself in propulsion system respons.^ [note Eq, (44)], The damping 
loss envelope components generally indicate that the major contributors to 
potential instabilities are associated with cmgine gimbal motion and oxidizer 
tank motion. 

Envelopes of total damping loss consisting of the sum of the com- 
ponents ^^tf ^^to presented in Figure 57 for various OMS fill 

fractions, indicr'^ing the potential for serious pogo destabilization over a 
wide range of model parameters. The initiation of roll off in damping loss 
with increasing tank modal amplitude depends on the OMS tank fill level. 

Near tank depletion, serious destab* lization is possible in tank dominated 
OMS pod modes > !)• The envelopes suggest that pogo suppression 

may be required for pod modes other than the localized engine gimbal mode. 

3,4 PRELIMINARY SUFPREoSOR STUDY 

The nominal and dispersed parameter stability results suggest that 
modification of both propellant circuits may be necessary due to the bi- 
propellant nature of the destabilizing system action (Section 3.3,3). Due to 
the potential for strong destabilization over large v \riation of pod struc- 
tural mode shapes (Section 3.3.5), a structural modification does not now 
appear to be as potentially effective as an accumulator modification in view 
of the large uncertainties present in the structural modes# 

3 3 

A suppressor configuration consisting of 82 ckx (5 in. ) compliant 
accumulators (i.e., negligible incrtance) at the fuel and oxidizer engine 
interfaces [corresponding to the contractual interface between Aerojet Linuid 
Rocket Company (ALRC) and McDonnell Douglas Astronautics Company 
(MDAC)] was chosen for preliminary evaluation. Fuel and oxiaizer circuit 
fundamental modes reduced to 22. 9 and 19.4 Hz, respectively, when sup- 
pressors were incorporated in the math model while the higher mode frequen- 
cies remained above 90 Hz. The calculated effect on the propulsion system 
dynamics was favorable above 20 Hz, (Chamber pressure response functions, 
altered by incorporation of accumulators, are illustrated in Figures 58 and 59 
along with the corresponding response functions without accumulators. 
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Figure 58. Effect of 82 cm (5 in. ) Accumulator on 
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Basically, the engine-excited response is phase stable above 20 Hz in spite 
of a general increase ir. magnitude, and the oxidizer tank-excited response is 
well attenuated above 20 Hz. The stability of the nominal primary axial 
mode is slightly degraded below JO Hz (minimum 8.9 dB stability margin) 
and enhanced above 20 Hz, and the engine gimbal mode at 45 Hz is phase- 
stabilized due to the suppressors. 

Damping loss envelopes have not been calculated for the OMS with 
accumulators. It is possible that significant destabilization may occur 
below 20 Hz. Evaluation of larger compliant accumulators, resistive ac- 
cumulators, and tuned accumulators should be pursued in further studies. 

3 

It is currently estimated that 330 cm (20 in.' ) compliant accumulators 
will reduce the fundamental propulsion system modes to about 10 Hz and 
enhance stability in the 10-20 Hz range. 

3.5 SUMMARY 

A preliminary mathematical model describing OMS intra-pod 
structural /propulsion system dynamics has been developed and analyzed 
for evaluation of pogo destabilization tendencies. A strong destabilization 
occurred in the localized 45-Hz axial engine gimbal mode that persisted at 
most OMS tank fill levels and under wide variation of the gimbal mode fre- 
quency. For nominal parameters in tb»» preliminary model, an instability 
occurred in this mode. Other potential destabilization tendencies were 
identified upon variation of pod structural mode shapes. Worst-case damping 
loss envelopes reveal the potential for serious destabilization in tank coupled 
modes over a wide range of tank interaction. 

3.3 

An initial investigation of small 82 cm (5 in. ) compliant accumu- 
lators located at the engine interface (i.e., the MDAC/ALRC interface) 
shows their benefit in suppressing the engine gimbal mode instability in the 
preliminary model. Degradation of stability at low frequencies (f < 20 Hz) is 

possible with such small accumulators an«l it may be necessary to increase 

3 3 

their volume to about 330 cm (20 in. ). 
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The present work represents the first OMS pogo stabilit r invest!-* 
gation« Future studies should include evaluation of (1) other feedsystem 
configurations (pod crossfeed, PBK supply, RCS supply), (2) lateral struc- 
ture/propulsion dynamics, (3) more extensive parameter sensitivity analy- 
sis, and (4) resistive and tuned accumulators. In addition^ pertinent OMS 
structure and propulsion system testing should be defined for identification 
of system dynamic parameters. 
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4. ERRORS IN FREQUENCY RESPONSE MEASUREMENTS 


Various hydraulic pulsing programs will be conducted to determine 
experimental frequency responses for identification of the linearized dyna- 
mics of propulsion system elements. These include the CTL-V program 
(LPOTP, interpump duct, accumulator), single engine tests, and the Main 
Propulsion Test (Sec. 4 of Ref. 11). 

The measurement of frequency response for these tests is difficult. 

If the induced system responses are too low, the results will be inaccurate 
due to the high self noise of the test system from operation of turbomachines/ 
thrust chambers, etc. If the induced system responses are too large, the 
results will be distorted by nonlinearities. In this section we estimate the 
random error in frequency response measurements due to system noise, and 
we consider other practical matters of significance to the accuracy of such 
measurements. 

A schematic diagram for the problem at hand is shown in Figure 60. 
It depicts the use of an external excitation i(t) to produce the two system 
responses x(t) and y(t). As shown these responses are made up of system 
random self-noise contributions m(t) and n(t), in combination with the re- 
sponse signals u(t) and v(t) which are puraly the result of the excitation i{t). 



Figure 60. Schematic of Measurement Situation 



We will be concerned with the use of sinusoidal excitation as a means to 
determine various frequency responses of the system. Reference is made 
frequently to Appendix V which was prepared by J.S. Bendat, a consultant for 
this investigation, 

RANDOM ERROR IN DETERMINATION OF H, 


4. 1 


2S}L 


As given in Eq. (F-13), the desired frequency response function 


H (f) can be estimated without bias by 
xy 


H (f) = G. (f)/c. (f) 
xy ly V IX 


(46) 


where and are cross- spectral density functions for the sub- 

scripted signals. We define the normalized mean-square error for the 
A 

magnitude of H (f) 




(47) 


Assuming Gaussian statistics, the fractional error in |H (f)| will be less 

xy 

than €(f) with 68 percent probability, and less than 2€(f) with 'IB percent 
probability. For the case of small errors and uncorrelated ..oise signals, 
the result is 


€^(f) « a^(f) + a^(f) (€^«1) (48) 

X y 

and in the worst case with correlated noise signals^ the result is 

€^(f) « [a^(f) 4 a^(f)]^ (€^ « 1) (49) 

X y 

where a^(f) = G (f)/u^{f)T(f) and a^(f) = G (f) /v^(f)T(f). The sinusoidal 
X mm / o y nn / o 

signal amplitudes u^(f) and v^(f) are at the frequency f. The analysis time for 
the sinusoidal signal at the frequency f is T(f). This result is comparable to 
Eq, (F-50) except in two respects; 
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a. The two signal/noise ratios are not equated here. 

b. The estimate of E[ |g. |^] in Eq. (F-I6) is halved based on 
the heuristic argumenf that, since i(t) is deterministic, the 
mean-square error in |G. | corresponds to that of the mean 
value of narrow-band noise. 


Note that the rms error t is independent of the bandwidth of the 
spectral analysis. An estimate of the error, however, is implicitly depen- 
dent on because of uncertainty in the estimates of the noise spectral 
densities G and G . Such uncertainty decreases with the B T. . , pro- 
duct [equivalently the number of average n^; see Eq. {F-47)], where 
is the total analysis time for the noise signals. This uncertainty is expected 
to be relatively small for the various pu’ jing test programs since 
should be relatively large. Fox example, a practical case is B^ = 0.5 Hz 

and T. , = 50 s, yielding a 20 percent uncertainty in G and G and a 
total « 7 e> r 7 

consequent 14 percent uncertainty in €. 

Another implicit influence of the bandwidth B is the fact that the 
A e 

estimate H^y,(f) will be some average over the bandwidth B^. Thus, the 

expected rate of change of H (f) with frequency should be considered in 

xy 

establishing the analysis bandwidth so as to avoid unacceptable smoothing. 

In the case of the measurement of the elements of the LPOTP transfer matrix 
(CTL-V testing), the resulting functions are quite slowly varying with fre- 
quency (see Figures 32 through 39) and a relatively wide analysis bandwidth 
is acceptable from this standpoint. 


4. 2 SWEEP TESTING 

One way to perform a frequency response test is to continuously 
change or sweep the excitation through the frequency range of interest. 
Another approach is to perform the testing at discrete values of frequency, 
holding the signal for some specified duration at each frequency. The time 
T(f) in the error expression of Eq. (48) i« the time that the excitation is 
maintained at the frequency f. For the svept case, T(f) is taken to be the 
time that the excitation lies within the an* lysis bandwidth B . 
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For a linear sweep (that is, a constant rate of change of frequency 
f©), we have 


T(0 = 


(50) 


The nunnber of cycles of excitation, Ng(f), occurring within the analysis 

bandwidth B (f) 
e 

Ng(f) = fT(f) = (51) 

showing that, for constant bandwidth and sweep rate, the nunnber of cycles is 
proportional to the frequency. 

For a logarithmic sweep (that is, rate of frequency change pro- 
portional to the frequency), we have f in Hz/s given by 


f = (log 2/60)R , f = 0.0116 R, f 
®e log log 


(52) 


where i® tl'® common expression of sweep rate in octaves per minute. 

Therefore, the analysis time in seconds is 


T(‘ 



= 86.6 B ())/r, { (B^/3f^«l) 

e / log e 


(53) 


where f is the mid -frequency of the band 3^(f). Moreover, the number of 
cycles of excitation within a bandwidth B^(fj is 


Nb«> = ®e'«/®log 


(B^/3f^<K 1) 
e 


(54) 


The frequency being varied linearly is the true instantaneous frequency, 
which is proportional to the time rate^ of change of excitation phase. Thus, 
the linear frequency variation is f^+f^t, the excitation must have the form 
sin 2rr[f^, + (fjj/2)t]t. 


if 
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Observe that, for a constant bandwidth, Ng is independent of frequency. 

For example, if = 1 Hz and ~ ^ oct/min: at f=2 Hz, the analysis 

time T is 14.4 sec; at f=40 Hz, T = 0. 72 sec; in both cases the number of 
cycles Ng = 29. 

4. 3 ACCURACY OF SINUSOIDAL AMPLITUDE 

It is clear that some minimum number of cycles Ng is required to 
obtain an accurate estimate of the sinusoidal amplitude of a signal. Use of 
straight integration to obtain the mean*square value of a sinusoid yields an 
upper bound on the normalized error of 1/47TN, where N is the number of 
cycles analyzed. To assure a fractional error of no more than in the esti- 
mate of the ratio of the amplitudes of two sinusoidal signals, the require- 
ment on N is 


N > 1/47T6 
o 


(55) 


For example, for at most a 2 percent error in the estimated amplitude 
ratio, N > 4 cycles. Such a requirement may not prove to be very con- 
straining for practical testing, but it should be kept in mind. 

4. 4 PLANNING A TEST 


Given estimates of the noise spectra (G and G ) based upon 

^ mm nn ^ 

prior noise measurements and a desired maximum rms error € for the 

m 

determination of set the conditions [based upon Eq. (48) assuming 

uncorrelated noise] 


a^(f) and a^(f) < /z 
x y m/ 


(56) 


There conditions translate into conditions upon the sinusoidal amplitudes 
and associated analysis time 

“>T(£) > ,57) 
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An example is now creattid to help crystallize the various error 

matters related to the planning of a test for accurate frequency response 

measurements. One system response x(t) is considered for this example. 

.2 

Suppose that x(t) is a pressure variable, and that = 1 psi /Hz and 

= 1 Hz over a 2 to 40 Hz frequency range of interest. Suppose also that 

€ =0.05 (maiximum rms error is 5 percent). The condition in Eq. (57) 

m 

yields 


u^(f)T(f) ^ 2(l)/(0.05)^ s 800 

Table 6 shows the combinations of amplitude and associated minimum analysis 
time in the first two columns. The next three columns show for the lowest 
frequency of 2Hz: Ng, the minimum number of cycles required within the 
bandwidth [from Ng(f) = fT(f) in Eq. (51)]; f^, the maximum constant sweep 
rate [from Eq. (50)]; and maximum logarithmic sweep rate [from 

Eq. (53)]. The iigiii-most three columns contain the same information for 
the highest frequency of interest of 40 Hz. 

Observe that a high-frequency, low-amplitude signal imposes a 
severe requirement for slowness of a logarithmic sweep. Note also that the 
low-frequency, high-amplitude condition yields a requirement for only four 
cycles v’ithin the bandwidth, which is apj caching the limit of concern for the 
accuracy of the sinusoidal amplitude itself [see Eq. (55)]. 

In planning an actual test, the various limits should be calculated 
based upon the best available pretest information for all of the system 
variables to be used for frequency response determinations and over the 
entire frequency range of interest. This should be done for several values 
of rms error. The results can then be employed to develop a lo^ cal test 
plan, based upon reasoned engineering judgment of all of the factors in- 
fluencing accuracy. 
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Table 6 


. Results for a Test Planning Example {€ * 0.05, 

B = 1 Hz, G =1 psi^/Hz) 
e mm 



f = 2 Hz 

f » 40 Hz 


Amplitude 

u psi 
o ^ 

Analysis Time 
T sec (>) 

cycles 

{>) 

• 

f 

o 

Hz/sec 

(<) 

log 

oct/min 

(<) 

^B 

eye le 8 
(>) 

• 

f 

o 

Hz/sec 

(<) 

Ri 

log 

oet/min 

(<) 

20 

2 

4 

0. 5 

21. 7 

80 

0. 5 

1. i 

10 

8 

16 

0. 13 

5.4 

320 

0. 13 

0.27 

5 

32 

64 

0.03 

1.4 

1280 

0. 03 

0.07 












4. 5 CONCLUSION AND RECOMMENDATION 

An engineering basis has been laid for the advance planning and 
evaluation of propulsion tests for frequency response determinations. It is 
recommended that the results developed herein be applied to the planning of 
the applicable Shuttle propulsion tests. 
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APPENDIX A 


EQUATIONS 1 MPS STABIUTY ANALYSIS 

A. 1 FLUID DYNA^lIC EQUATIONS FOR 
PROPULSION SYSTEM 

The fluid dynamic equations for the elements of the propulsion sys- 
tem model are given below. 

First Feedline Segment: 


p 

2 " “ 

(A-1) 

(ItsZjCjo 

'^2 - '’4 - ^2°t = » 

(A-2) 

Downcomer Segment: 



^6 ■ “l C4 ■ °I2 

ht - ^2 - *':24''4 1 ' 0 

(A-3) 

- “2P4 - “22 

|q,-»C,„P2-.C24P4| = 0 

(A-4) 


owncomer Corner to Manifold: 


(l+sZ^Cjj)P^ - Pg - Z^Qg = 0 (A-5) 

Ut.Z5Cj2)F8-P|o-/-5|c>5-=C,,pJ =0 (/.-6) 

(I + |q^ - - sCj^Pgj *0 (A-7) 

Manifold Continuity: 


A-1 



Manifold to Upper Engine LPOTP: 


(l+sZ^Cj^) - Pj 4 - Z^Qj 4 = 0 


(A- 9) 

(l*»ZgC, 5 ) - P,g - Zg|Q,^ - sC,gP,g 

1 * ^ 

(A. 10) 

*"l 6 ■ ^17 ■ ^9l®14 ■ ■ “''15^.41 

= 0 

(A- 11) 

Manifold to Lower Engine LPOTP: 

(1 +sZjqCj^) Pj 2 - Pjg - ZjgQjg = 0 


(A- 12) 

(1 +sZjjCj^) Pjg - P,Q - Zj J Qjg - sCjgPj 2 

1 = 0 

(A-13) 

^20 ■ ^21 ■ ^12 1^18 ■ ®^16^12 ■ ®^17^isl 

S 0 

(A- 14) 

Upper Engine LPOTP: 

^40 ■*■ ®^pl^I7 ■ ^14 ^ ® 1^3^40 ■'■ ^14^12 ^15^1 

(1 

0 

{A- 15) 

^40 ■ ("'i + lJPn + ^/p^40 ~ ® 


(A- 16) 

Upper Engine Inter-Pump Duct: 

( 1 +sZj 3 Cjg)P 4 Q - P^j - = 0 


(A- 17) 

{l + sZj4Cjg)P^j - P^2 - Z 41 IQ 40 " 18^40 1 

= 0 

(A-18) 

(1 +sZj5C2o)P42 - P 43 - ^ 42(^40 ■ 18^40 ‘ *‘"19^41^ “ 

0 

(A- 19) 

'43^a^^43 ■ ^44 ‘ ^431^40 ‘ ®^18^40 “ 19^41 " 

“'^20^42 

1 = 0 


(A-20) 
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Upper Engine HPOTP: 


^44 ®^p2^44 ■ ^40 ^a^43 ' ® ^18^40 ■*■ ^19^41 ■•■ *^20^42 


0 

(A-21) 


^45 " ^*^2 ■'■ ^^^44 ^ ^hp^44 “ ° 


(A-22) 


Upper Engine Injector and Chamber: 


P 45 - l«c"V°44 = » 


1> , - R Q. . = 0 
cl c 44 


(A-23) 

(A-24) 


Lower Engine LPOTP: 

^50 ®^pl^21 ■ ®18 ^ ® 1^3^50 ^16^12 ^17^18 I (A-25) 


^50 " ^"^"^^^^^21 ^/ p ^50 ■ ® 


(A-26) 


Lower Engine Inter-Pump Duct: 


^'"®2i7S.i)P50-P51-^50Q50 = ° 


(1 +sZjgC 23 ) P52 - P53 - 252,053 - SC 2 JP 5 Q - SO22P5J = 0 


(A- 27) 


(1 +sZjgC 22 )P 5 i - P 52 - Z 5 Jq^q - sC 2 jP 5 ,^ *0 (A-28) 


(A-29) 


^53^a| ^53 ' ^54 ‘ ^53 PsO " ®^21^50 ’ ®^22^51 “ ®^23^52 " ® 


(A-30) 
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Upper Engine HPOTP: 


^54 ®^p2^54 ■ ^50 ^ ^a^53 ^ ® 1^21^50 ^ ^22^51 ^ ^23^52 | ^ ° 

(A-31) 

P 55 -(m 2 + l)P 54 .Z^pQ 34 = 0 (A-32) 

Lower Engine Injector and Chamber: 

P 55 - (R,. ♦ ZjlQj^ = 0 (A-33) 

Pc2 - «c«54 = “ 

where the P and Q denote the pressures and absolute flows at different 
m n ^ 

stations in the propulsion system (see Figure A-1). The coefficients 

denote the local compliance values that are employed to represent the effect 

of compressibility in related elements of the feedline; the locations of these 

compliances are shown in Figure A -2. The coefficients Z denote the 

m 

impedance associated with the various sections of the feedline. This im- 
pedance is defined by 


Z — Ls + R 


where L^ R denote the inertance and resistance associated with the line seg- 
ment. The terms ^p 2 * ^*^1 ^^2 ^ pump inlet cavi* 

tation compliance and p\imp gain associated with the LPOTP and HPOTP, 
respectively; the coefficient denotes the resistance of the thrust chamber. 
The characteristics of the accumulators are defined by 


Y 

a 


= (l^) +6.^ j 


(A-35) 


A-4 




where and denote the inertance, frequency (rad/sec), and damping 

of the accumulator, respectively. The frequency and damping ratio are 
defined by 

^a = (A- 36) 


where C and R denote the compliance and resistance of the accumulator. 

fll 

The coefficients a.. are associated with the compressible flow solutions for 
the downcomer and are defined as follows: 

“1 1 ~ “22 ■ ^ 

(A-37) 

^ sinh B 

12 * ■ ^3 B 

(A-38) 

„ sinhd 

Zj 

(A-39) 

where the parameter B is defined by 

/ R 

(A-40) 


where i is the travel time for an acoustic wave in the downcomer segment; 

L 3 and denote the inertance and resistance associated with the downcomer. 

A 

The motion of the structure X is represented by a series of normal 

modes 


X $„(r) 


{A-41) 


n 


th ^ 

■ here q (t) is the n generalized coordinate and 9 (r) is the associated mode 

“ th ^ 

shape. The equation of motion for the k generalized coordinate is 

= ^^(tb)Q^ 4 P,A,«W(tb) 


(A-42) 


A -6 



where (tb) ia the modal tank-bottom pressure and Q is the relative 
volume outflow from the propellant tank. The tank bottom pressure, P^, 
is related to the vehicle motion by 

n 

The F. comprises the drag forces on the feedline segments, the inter-pump 
ducts and discharge lines, the forces at the feedline corners, the forces on 
the LPOTP and HPOTP, and the forces on the thrtist chambers. 
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APPENDIX B 


NUMERICAL DATA FOR MPS STABIUTY ANALYSIS 


PROPUL5ION SYSTEM 

PARAMETERS 


Resistance, 

MN s/m^ 

(10"^ sec/in, 


0.022 

( 1.32) 

■<3 

0.0<)7 

( 3.97) 


0 . 078 

{ 4.62) 


0.029 

( 1.68) 

R? 

0.031 

( 1.81) 

^8 

0.027 

( 1.61) 

*^10 

0.031 

{ 1.81) 

«11 

0.027 

{ 1.61) 

^40 

0.40 

(23.5 ) 

«41 

0.78 

(46.0 ) 

^42 

0.41 

(24.0 ) 

R 43 

0.016 

( 9.6 ) 

^50 

0.40 

(23.5 ) 

«51 

0.78 

(46.0 ) 

^52 

0.41 

(24.0 ) 

R 53 

0.016 

( 9.6 ) 

R 

c 

42.3 

(2500) 

"/p 

6.8 

( 400) 

"hp 

44. 7 

(2640) 

R. 

J 

68.6 

(4050) 


Resistance in SI units is based upon pressure divided by volume flow; in 
English units, weight flow is employed. 
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Inertance, 


(10~^ sec^/in.^) 


10 “ ^ MN 3^/m^ 


^1 

3. 7 

(0.22) 

^2 

32. 2 

(1.9 ) 


220. 

(13) 


22. 0 

(1.3 ) 

^5 

14. 4 

(0.85) 


5. 1 

(0.3 ) 

^7 

22. 0 

(1.3 ) 

^8 

22. 0 

(1.3 ) 

^10 

23. 7 

(1.4 ) 

hi 

20.4 

(1.2 ) 

ho 

104. 

'^.0 ) 

hi 

o2. 7 

(3.7 ) 

hz 

i3. 9 

(2.0 ) 

hs 

28. 8 

(1.7 ) 

ho 

104. 

(6.0 ) 

hi 

02. 7 

(3.7 ) 

hz 

33. 9 

(2.0 ) 

hs 

00 

• 

00 

(1.7 ) 

hp 

17. 0 

(1.0 ) 

hp 

42. 4 

(2.5 ) 

L. 

J 

84, 7 

(5.0 ) 

h 

18. 6 

(1.1 ) 


Inertance in SI units is based upon pressure divided by volume flow; in 
English units, weight flow is employed. 
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Compliance, 

lO'^^m^/MN 

o 

1 

• 

""lO 


9.6 

(16.0 ) 

^11 


5. 6 

( 9.4 ) 

^12 


9.6 

(16.0 ) 

^13 


5.Z 

( 8.6 ) 

^14 


1,8 

( 3.0 ) 

^15 


2.3 

( 3.8 ) 



1.9 

( 3.2 ) 

^17 


2. 3 

( 3.8 ) 

^18 


0.36 

( 0.6 ) 

^19 


0.6 

( 1.0 ) 

"^20 


0.45 

( 0.75) 

Si 


0. 36 

( 0.6 ) 

^22 


0.6 

( 1.0 ) 

^23 


C 45 

( 0.75) 

^24 


9.6 

(16.0 ) 

c , 1 

pi 

See Table 

1 


c , 

p2 





Pump Gain 


m j + 1 
m^+l 


See Table 1 


Transit Time, sec 
T = 0.0657 


Compliance in SI units is a volume change per unit pressure change; in 
English units, it is the weight of liquid displaced per unit pressure change. 
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1 


1 


I 


I 


Areas 

in-2 2 

10 m 


^1 

14. 7 

227 

^2 

7.2 

112 

Aj 

2.0 

31 

^4 

0.81 

12.6 

^5 

10. 2 

158 


B.2 STRUCTURAL MODE DATA 

Modal data are provided in the table for the second (L^) and fifty- 
fourth modes at liftoff. These inodes proved to be of most interest in 

the stability calculations. 
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1 1 


a 


The modal amplitudes in the table are normalized to a generalized 
4 2 

mass of 1. 1453 X 10 lb-sec /in. In the table, (tb) refers to the modal 

2 ^ 3 

tank-bottom pressure (givren m units of lb-sec /in. ); 0(e^) and ^(e^) denote 
the modal amplitudes at the upper and lower engine gimbal blocks, respec- 
tively. 
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APPENDIX C 

SIMPUFIED OMS STRUCTURA L MODEL DETAILS 


Included in this appendix are the 11 -DOF OMS mass and stiffness 
matrices. Tables C-1 and C-2, which are primarily based on RI/SD data, 

A summary of modal parametc^rs for the primary axial modes and the local 
engine gimbal mode is presented in Table C-3, 
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Table C-1. OMS Pod Masses 


DOF 

Description 

1 

RCS fuel tank - X 

2 

RCS ox tank - X 

3 

OMS fuel tank - X 

4 

OMS ox tank - X 

5 

OMS H tank - X 
e 

6 

7 

PCS thrusters eg - X 
Engine X 

8 

Engine Y 

9 

Engine Z 

10 

Engine 

11 

Engine 9 
^ z 


fi denotes OMS fractional fill level 
V denotes RCS fractional fill level 
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Table C-3, OMS Structural Mode Summary 


Fill Condition (%) 

Primary Axial Mode 

Engine Gimbal Mode 

OMS 

RCS 

Hz 

M (Ibm) 


^to 

— 

Hz 

r 

M (Ibm) 

‘Oti 


100 

100 

10. 35 

5. 82x10"^ 


2. 155 


284 

-0. 009 

-0.003 


50 

10.51 

6. 10x10^ 

1. 708 

2. 312 


284 

-0. 008 

-0.003 


1 

10.60 

6.42x10^ 

1. 716 

2.430 


372 

-0. 007 

-0.001 

75 

100 

11.57 

4. 14x10^ 

1. 631 

1. 922 


284 

-0. 012 

-0. 004 


50 

1 1. 86 

4. 25x10^ 

1 . 657 

2. 137 

44. 87 

284 

- 0 . on 

-0. 004 


1 

12.02 

4.51x10^ 

1.669 

2. 288 

45. 03 

357 

-0.009 

-0.001 

50 

100 

13. 05 

3. 05x10"* 

1. 509 

1. 533 

44. 87 

284 

-0. 017 

-0. 006 


50 

n. 86 

2. 54x10^ 

1. 562 

1. 836 

44. 88 

283 

-0. 017 

-0.005 


1 

14.20 

3.22x10"* 

1. 579 

2. 037 

45. 04 

352 

-0. 013 

-0. 002 

! 25 

100 

17.62 

2. 21x10"* 

1.433 

1. 139 

44. 93 

283 

-0.034 

-0.010 

1 

I 

50 

16. 85 

1. 28x10"* 

1. 340 

1. 267 

44. 94 

282 

-0. 033 

-0. 010 


1 

18.08 

1.06x10'* 

1.359 

1. 500 

45. 07 

338 

-0. 027 

-0. 003 

10 

100 

25. 10 

6. 55x10^ 

0.769 

1. 120 

45. 10 

282. 

-0. 081 

-0.018 


50 

22. 50 

1. 06x10^ 

1. 066 

0.684 

45. 10 

281 

-0. 079 

-0.017 


1 

21.87 

3.88x10^ 

1. 026 

0. 874 

45. 16 

282 


-0.007 

3 

100 

26. 16 

4.49x10^ 

0. 547 

0.472 

45. 58 

303 

mgm 

+0.002 


50 

23.09 

7. 84x10^ 

0. 884 

0.457 

45. 58 

342 


+0. 002 


1 

23. 54 

3. 17x103 

0. 807 

0. 570 

45. 55 

296 

HI 

+0.003 

















APPENDIX D 


OMS PROI»ULSION MODEL DETAILS 


D. 1 FEEDUNES 

In general, the allocation of nodal inertance, resistance and stiffness 
is accomplished in a straightforward manner based on line geometry, pres- 
sure budget and reduced bulk modulus data (accounting for pipe radial elas- 
ticity) presented in Tables D- I and D-2. Specific complications associated 
with free bellows and engine manifold geometry, however, required special 
consideration . 

The engine fuel line segment JKL and oxiJizer line segment KL, 
respectively, have free bellovs sections restrained by adjacent curved pipe 
sections cantilevered at the engine valve housing. The increased volumetric 
compliance due to the bellows is approximated on the basis of flexibility of 
the adjacent pipe sections illustrate<l in Figure D-1. The complementary 
strain energy due to bending is 


e- 



o 



2 3 2 

A^R^P^ ,, ^,2 

— (1-cos^) dO 


(D-1) 


with 0“^ equal to jr/2 and 7T for the fuel and oxidizer circuits, respectively. In 
accordance with Castigliano's theorem, the volumetric compliance due to 
bellows is 



belloA^s 


1 


Su 

c 




(D-2) 


D. 1 



Table D-1. OMS Fuel Circuit Hydraulic Properties 


Line = 8. 138 X lO"^ Ib-sec^/in.^) 


Segment 

O.D. 

(in.) 

t (in.) 

L (in.) 

A (in,^) 

5 

Bulk Modulus (10 pai) 

Acceleration 
Head (in.) 

/ip (psi) 

AB 

1. 

5 

0. 03 

11. 3 

1 . 6286 

1. 175 

- 1 1. 3 (tank) 


BC 





20. 0 

3.2572* 



. 

4. 0 

CD 





14. 0 

1.6286 



-14. 0 (tank) 

M M M M 

DE 





4. 0 





- 4. 0 (tank) 

M M — «■ 

EF 





25. 0 





.... 


FG 





28.0 





-28.0 (gimbal) 


GH 





26. 15 





- - - - 

24. 0 

HI 





13. 2 





+ 9. 6 (gimbal) 

M M 

IJ 

1 

5 

0. 03 

8. 37 

1.6286 



.... 

4.0 

JKL 

1 

25 

0.025 

12. 79 

1. 131 

1. 175 

-10. 9 (gimbal) 

m m m mm 

DE' 

2 


0. 035 

7, 4 

2. 9255 

X. 158 

.... 

- - - - 

E'F' 

2 


0. 03' 

25. 0 

2. 9255 

1. 158 

+25. 0 (tank) 


F’G’ 

1 

5 

~0. 03 

10. 0 

3.2572* 

1. 175 

• MM 



2, Engine Manifold/Regen/Injector 


Point/ 

Segment 

2 5 

Inertance (lb-sec /in, ) 

3 

Volume (in. ) 

5 

Bulk Modulus (10 psi) 

Acceleration 
Head (in. ) 

dP (psi) 

M/MN 

2. 85 X 10’"^ 

37.0 

1. 5 

-4.0 (gimbal) 

M M M M 

N/NO 

9.75 X 10-4 

37.0 

1. 5 

+ 11. 3 (gimbal) 

12.5 

O/OP 

1. 06 X 10-3 

6.5** 

1. 5 

+ 11.3 (gimbal) 

73. 5 

P 

2. 70 X 10-4 

N.A. 

N.A. 

“ 

7.0 


Tv/o lines in parallel at double valve 

'^'^Estimate of stiffness includes injector structure (AlP/Aw*^\0^ psi/in*^) 





























Table D-2. 


OMS Oxidizer Circuit Hydraulic Properties 


1 . 


Line (p^ = 1.3525 X 


lO”^ Ib-sec^/in.^) 


Segment 

O. D. 

^—1 
(in. ) 

t (in. ) 

t (in. ) 

A (in.^) 

Bulk Modulus (10^ psi) 

Acceleration 
Head (in.) 

AP (psi) 

AB 

1. 

5 

0, 03 

9.6 

1.6286 

1. 060 

- 9s 6 (tank) 

— 

1 

1 

1 

BC 





20. 0 

3.2572* 




6, 0 

CD 





19.6 

1.6286 



-16.0 (tank) 

« - . • 

DE 





4.0 





- 4, 0 (tank) 

- - - - 

EF 





5.4 





... - 

.... 

FG 





28. 0 





-28.0 (gimbal) 

-- -- 

GH 





29. 445 





m m ^ 


HI 





13. 1 





+ 13. 1 (gimbal) 

---- 

IJ 

1. 

5 

0, 03 

11. 035 

1.6286 



- - - - 

7.0 

JK 

1. 

25 

0. 025 

10. 9 

1. 131 



+ 7. 02 (gimbal) 

m ^ ^ m 

KL 

1. 

25 

0.025 

11. 13 

1. 131 

1. 060 

- 2, 83 (gimbal) 


DE' 

2. 


0.035 

6. 0 

2. 9255 

1. 047 

- - - - 


E'F' 

2, 


0.035 

20.0 

2. 9255 

1. 047 

+20. (tank) 


F'G' 

1. 

5 

~0.03 

10. 0 

3.2572* 

1.060 


— • — 


2. Engine Manifold/Injector 


Point/ 

2 5 

3 

5 

Acceleration 

AP (psi) 

Segment 

Inertance (lb- sec /in. ) 

Volume (in. ) 

Bulk Modulus (10 psi) 

Head (in. ) 

M/MN 

3. 96 X lO'"^ 

6.5** 

1. 3 

» « - « 

71.0 

N 

4.49 X 10-4 

N. A. 

N. A. 

• - — 

7.0 


Two lines in parallel at double valve 

4 3 

Estimrte of stiffness includes injector structure (LPlAV 10 psi/in. ) 






















Figure D-1. Free Bellows Dilatational Compliftiic'j Models 
^ (R s 3. 5 in., t s 0. 025 in., d a 1. 2s in. , 

E s 3 X 10^ pal) 




(a) FUEL 


K 

(b‘ OXIDIZER 



For the fuel and oxidizer circuits, the respective values are 



= 3.639 X 10"^ in.^/psi 


4,815 X 10”^ in.^/psi 


(D-3) 


Thus, the effective hydraulic stiffnesses of the fuel line section JKL and 
oxidizer line section KL are aj>proximated by 


K 


JKL 


fuel 




OX 


r 1 A 

*FKL JKL (9 



® JKL ^ 

fuel bellows i 

6,269 X 10^ psi/in.^ 



r / A 

''KI- KL , / 9 \ 

1 

-1 

B 

1 KL 

ox bellows 



= 1. 663 X 10^ psi/in.^ 


(D-4) 


(D-5) 


While the bellows produce significant local stiffness reduction in the feed 
system stiffness matrices, the effect on the propulsion system node shapes 
and frequencies is small. This is due to the fact that the bel.ows are close 
to the chamber interfaces. 

D. 2 ENGINE MANIFOLD MODEL 

The engine manifold assembly illustrated in Figure D-2 consists of a 
short fuel entry duct, a lower fuel distribution manifold, a fuel regenerative 
cooling assembly, upper fuel and oxidizer distribution manifolds, and fuel 
and oxidizer injector faces. A preliminary examination of the structural wall 
thicknesses and overall dimensions intUcates that the dilatational structural 
compliance of the engine manifold is negligible. Therefore, compliance of 
the engine propellant circuitry is mainly due to the fluid bulk moduli. 


D-6 




^26.4 cm_^ 
(10.4 in.) 



42.2 cm - 
(16.6 In.) 


R s: 18cm (7 in.) 

Agyg • 5.4cm^ (0.83 In. ^) 


Figure D-2. OME Manifold Schematic 


D-(. 




Calculation of the resi'^tance and stiffness coefficients for the engine 
manifold model was based on /olumt; and steady pressure ''•rep data indicated 
in Tables D-1 and D-2. In contrast, estimation of the engine manifold iner- 
tance distribution was complicated by geometric details m spite of the general 
assumption of one-dimensional flow. A summary of compone-'ii inertance 
values and individual methods of calculation are presented in Table D-3. The 
one component requiring spec al analytical treatment is the lower fuel distri- 
bution manifold. 


On the assumption of uniform leakage flow into the regenerative cooling 

chamber, the circumferential manifold flow distribution is related to the fuel 

entry duct flow, &. , as 
^ in 


GiO) 




(D-6) 


The perturbation in manifold Ilow kinetic energy, assuming a cross-sectional 
area distribution given by 


is therefore 



(D-7) 


T 



7T 



A(d) 


A(9) Rdd 


(D-8) 


This results in the following expression for lower fuel manifold inertance 


L lower fuel 
manifold 




h0 


in 


7T PR 
■ 2 A 


avg 




(D-9) 


and numerical integration of tliis equat on incorporating geometric data 
presented in Figure D-2 results in the inertance estimate given in Table D-3. 





Table D-3. OME Manifold Inertance Distribution 


Component 

Inertance 

(Ib-sec^/in,^) 

Head (in. ) 

Method of Calculation 

Fuel Entry Duct 

2.7 X 10“"^ 

- 4.0 

one-dimensional flow, de- 
tailed dimensions from 
ALRC drawing 

Lower Fuel Manifold 

3 X lO"'^ 

0 

circumferential flow with 
uniform leakage to re- 
generative cooling channels, 
volume ^"36 in.^ 

Regenerative Cooling Chamber 

1.65 X 10"^ 

22.4 

one-dimensional flow uni- 
form distribution among 
channels, detailed dimen- 
sions from ALRC data, 
volume ^31 in.^ 

Regenerative Chamber to Upper Fuel 
Manifold Interface 

1.33 X 10’^ 

0 

119 holes, length ~ 0. 3 in,, 
diameter ~ 0, 14 ir., uni- 
form flow distribution 

Upper Fuel Manifold 

4. 5 X lO”"^ 

0 

3 

volume -w6,5 in, , effective 
flow path length ^6 in. 

Fuel Injector Face 

4. 59 X 10"^ 

0 

based on ALRC supplied 
value for effective length/ 
area 


Oxidizer circuit components differ by P^l present analysis. 




D. 3 COMBUSTION DYNAMICS - APPROXIMATE 
STATE EQUATIONS 


Since the present formulation and analysis approach requires a set of 
first-order, real, constant-coefficient differential equations, an approxima- 
tion of Eq. (28) has been derived. The approximation was first used by us 
during Delta Stage II pogo studies in 1975. Relatively long time delays with 
T^f greater than were originally anticipated: thus, Eq. (28) was first 
recast in the form 

— 7* S 

(1 + T s)P^ = e (D-10) 

r c c 


P = R 0. + e af ao R 0 (D-11) 

c CO lO cf xf 

The exponentials were then approximated in terms of fourth-order Fade ex- 
pansions (Ref. 20) as 


e 


-P 


1680 - 840p 180p^ - 20p^ + p^ 

1680 + 840p + 180p^ + 20p^ + p'^ 


(D-12) 


Substitution of the expansion into Eqs. (D-10) and (D-11) results in a pair of 
higher-order differential equations. Finally, incorporation of the state 
variable transformation described in Ref. 21 results in the set of first-order 
differential equations shown b<*low. The approximation is extremely accurate 
for the range of parameters of interest. 


D-0 



‘3 


-1680 -840 

-180 -20 

'^do ‘^dr’do 

^•ir'^do ‘^df"‘^do 

.40 


df do 


800 


’^df'^do 


-10480 


^df" 


99200 
■^dr '^do 



R /r 
CO r 

-40 R^f/lT^rW 
800 

-10480 R,f/(T^rW 

992 00 

-40 a /r, 

CO do 

800 R Jt. 
cf do 

800 R /r. 

CO do 

• 10480 R Jt. 

cf do 

-10480 P /r. 

CO do 

99200 U Jt. 

cf do 

99200 R Jt. 

CO do 


do 


-1680 >840 

^do 


do 


>180 


do 


do 


-20 

*^lo 


9 . 


if 


G. 


(D.I3) 


D. 4 CLOSED- LOOP STABILITY ANALYSIS CONSTRAINTS 

Details of the constraint matrix employed in assembly of the coupled 
propulsion/structure system model are presented below in Table D-4. 


D-IO 



D-ll 


0BI6INAL PAGE fb 
OP POOR QUALITY 


Analysis Constraints 


Constrained Veriabie 


Name 

J_ 

til 

®fi 

1 

1. 

0f2 

2 

1. 


3 

1, 

®£4 

4 

1. 

^£5 

5 

1. 

®£6 

6 

1. 

®£7 

7 

1. 

®£8 

8 

1. 

X 

e 

21 

^t£ 

X 

e 

21 

^£ 

X 

21 

1. 

e 

e . 

12 

1. 

Ol 


13 

1. 

®o3 

14 

1. 

®o4 

15 

1. 

®o5 

16 

1. 

®o6 

17 

1. 

e - 

18 


o ( 

®o8 

19 

1. 

X 

e 

21 

^to 

X 

e 

21 

^to 

X 

21 


e 

®£8 

8 

1. 

®£9 

9 

1. 

®£10 

10 

1. 

e... 

11 

1. 

til 

®o8 

19 

1. 


20 

1. 

X 

21 

1. 

e 

X 

21 

1. 

e 


21 

^t£ 

X 

c 

21 

^to 

^£1 

1 

1. 

®ol 

ll 

1. 
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APPENDIX E 


SOLUTION OF THE FIRST-ORDER DYNAMIC EQUATIONS 


The dynamics of the OMS structural/propulsion system are described 
in terms of a set of first-order differential equations with constant coeffi- 
cients of general form 


|x| - A |x| = B |f| (E-1) 

The matrix A is square and B is generally rectangular. For closed-loop 
stability analysis, the external forcing function j fJ is null, and for propul- 
sion system frequency response analysis, the forcing function c‘>nsists of 
structural interface accelerations. 

The general homogeneous solution of Eq. (E-1) 

|x I = j X I (E-2) 

i 

\ 

is obtained by QR iteration (Ref. 22) which calculates the complete set of 
complex eigenvalues and eigenvectors. The set of column eigenvectors (0) 
(generally complex), if numerically perfect, orthogonalize the dynamic equa- 
tions as 

(\) Ul " ° 

v/ith 

|x| = 0|^| (E-3b) 


i 


E-1 



Evaluation of the numerical accuracy of the complex modal solution is 
accomplished by application of the Gerschgorin theorem (Rtf. 22' which yields 
the following information: 

a. Diagonal entries in 0 ^A<j> represent improved estimates of the 
eigenvalues, X 

b. Off-diagonal row and column sums 


■•k ’ Yl, % 

j = l,n 
j«k 

'k ' Y 

i=l, n 
i^k 


(E-4a) 


(E-4b) 


yield a radius R, about 

X, must lie. ' 

k 


kk 


within which the true eivenvalue 


The maximum radius 


R. - max(r, ,c.) (E-5) 

K K K 

is used to assess accuracy of the eigenvalue, The overall accuracy o.' 

the numerical eigenvalue solutions is evaluated based on the smallness of the 
radii and the shifts 

Modal forced response of Eq. (E-1) is conveniently calculated with 
knowledge of the homogeneous solutions. In the present study, frequency 
response subject to the harmonic forcing function 

■ 1 1^0 I 



E-2 





I 


is required for study of propulsion system dynamics and open-loop stability* 
The steady state modal dynamic solution of Eq. (E-1) incorporating Eq. (E-3b) 
is 



(E-7a) 


Although truncated modal solutions may be calculated, the frequ ncy responses 
calculwtted in the present jtudy include the contribution of all n des. 


E-3 
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APPENDIX F 


ERROR ANALYSIS OF TRANSFER FUNCTION MEASUREMENTS* 


This is a summary of work conducted to analyze errors in transfer 
function (frequency response) measurements based upon noise/signal con- 
siderations and coherence functions between the measured data. Quantitative 
formulas are obtained that give the standard error in transfer function esti> 
mates obtained by using an external excitation signal. Proper interpretation 
of measured results depends upon evaluation of quantities discussed in this 
appendix. 

F. J GENERAL DIAGRAM AND NOTATION 

A diagram is shown in Figure F-1 of a general system \\diose transfer 
function H (f) is to be determined. Quantities in Figure F-1 have the 
following meanings; 

i(t) = measured external excitation signal 
u(t) = input signal caused by i(t) 
v(t) = output signal caused by i(t) 
m(t) = input self-noise independent of i(t) 
n(t) = output self-noise independent of i(t) 
x(t) = measured input signal = u(t) + m(t) 
y(t) = measured output signal = v(t) + n(t) 

A special case of Figure F-1 occurs when i(t) = 0, corresponding to self- 
noise only, where x(t) = m(t) and y(t) = n(t). For this case, the transfer func- 
tion H (f) will be denoted by II (f). 
xy ' mn 


The author of this appendix is J.S. Bendat, acting as a consultant to 
The Aerospace Corporation. 

An example of application to the CTL-V testing: i(t) is the electrical signal 
to the suction pulser; x(t) and y(t) are the inlet pressure and flow to the 
LPOTP, respectively. The transfer function is one of those required for 
four-terminal definition of the LPOTP. 


F-1 




Figure F-1, System with External Excitation 

From the measured signals^ i(t)» x(t) and y(t), autospectral (power 
spectral) density functions, cross -spectral density functions, and coherence 
functions can be computed. These will be denoted by: 

G^^(f) = autospectral density function of x(t) 

G (f) = cross-spectral density lunction between x(t) and y(t) 

xy 

(f) = coherence function between x(t) and y(t) 

xy 

with similar notations for G (f), G..{f), G. (f), G. (f), y\ (f) andy? (f). The 

yy 11 IX ^ ly ' ix ly 

coherence function between x(t) and y(t) is defined by 


(f) =|g (f)lvG (f)G (f) (F-] 

xy ■ xy ■ / XX yy 

Z Z 

with similar definitions for y. (f) and y. (f). 

IX ly 

2 

As is well known, for any frequency f, y (f) is bounded between 0 

xy 

and 1 and is a measure of an ideal (constant parameter) linear relationship 
between x(t) and y(t). 


F-Z 




The dependence upon f will now be omitted in all the following for- 
mulas for simplicity in notation, but should not be forgotten. For Fi«rure F-l, 
since the self-noise terms m(t) and n(t) are independent of u(t) and v(t) 

G =G +G , G sG +G , G =G +G 
XX uu mm yy vv nn' xy uv mn 

(F-2) 

Define the input noise -to-signal ratio o and the output noise-to- signal ratio 
/?by 


a = G /G , ij = G /G 
mm uu nn \ 


(F-3) 


Then 


G = G (1 + a), G = G (1+^ 

XX uu ’ yy vv ^ 


(^- 4 ) 


;he coherence function y is given bv 

xy ^ 


= I G +G iVg G (l + a)(l ffl) 
xy I uv mn| / uu vv ^ 


(F-5) 


Thus, the coherence function is clearly a function of a and as well as the 
other terms shown, 

F. 2 SINGLE INPUT /SINGLE OUTPUT MODE US 

For simple single-input/single-output models, where x{t) is a meas- 
ured stationary random input signal and y(t) is its associated measured output, 
the optimum constant parameter linear system between x(t) and y(t) has a 
transfer (frequency response! function given by 


H = G /G 
xy xy XX 


F- ? 


(F-6) 



In polar form 


H 


xy 




(F-7‘, 


with gain factor jH^^j and phase factor 

The output autospectral (power spectral) density function can be 
decomposed into two parts (see Figure F-2) 


where 


G 

yy 


H 

xy 


2 


G 


XX 


+ G 

zz 


(F-8) 


H 


xy 


G = G 

XX xy yy 


(F-9) 


is called the coherent output powt-r spectrum and represents the portion of 
y(t) due to x(t). The second term 


G = ( 1 - )g 
ZZ \ xy/ yy 


represents the portion of y(t) due to unknown z(t) caused by extraneous noise, 

nonlinear effects, etc. Note that G approaches zero as the coherence func- 
2 z z 

tion approaches one. 


i(t) 

I 

I 

I 

♦ 

xtt) 


z(t) 



- y(t) 


Figure F-2. Single-Input/Single-Output Model 


F -4 




Estimates of found by means of Equation (F-6' which do 

not involve the use of any external excitation signal i(t). These estimates 
will always be biased if extraneous noise is present at the input. By using i(t), 
however, it is possible to obtain unbiased estimates of This will now 

be shown* 

F. 3 TRANSFER FUNCTION ESTIMATES USING 
EXTERNAL EXCITATION 

The Fourier transform relation of Figure F-2, indicated by capital 
letters, is given by 


Y = H X + Z (F-11) 

xy 

where all terms in Equation (F-11) are functions of frequency. When the 
external excitation signal l(t) is availr^ble. Equation (F-11) yields the result 


G, s H G. + G. 
ly xy IX iz 


(F-12) 


Hence if the "expected value" of G^^ (denoted by E[G ^^] ) ^ 0, which is a 
reasonable expectation in these applications, then 



= G. /G. 
ly ix 


(F-13) 


will provide an unbiased estimate of ^^y* From Equation (F-12) observe that 

H = H 4 (g. /G. ) (F-14) 

xy xy \ iz IX / 

Equation (F-13) will be used instead of Equation (F-6) to obtain estimates of 

the desired transfer function H . The variance error in these estimates 

xy 

will now be indicated. 
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Assume that 


is essentially a large con- 


stant* Then 




where 



E 

1 A 1 2 

H -H 

= E 

|g. /g. 

1 

1 xy xy| 


1 IZ IX 1 



(F-15) 


It is also reasonable to assume that estimates of G. will be obtained by 

^ iz 

averaging over n^ independent estimates," It then follows that (Ref, 23, p, 185) 


E 



< 


G..G /n, 
ii zz d 


(F-16) 


The quantity 


|G. I ^ G..G 

IX| IX 11 XX 


Hence, Equation (F-15) becomes 


{I"xy-%0'S 


A2 

r 


where 


(F-17) 


(F-18) 


a2 

r 


G /n.y. G 
zz d IX XX 



)g /n 
/ yy/ 


,/g 

I IX XX 


(F-19) 


using Equation (F-10) for G . Note that all terms in Equation (F-19) are 

z z 

measurable quantities. 

Geometrically, Equation (F-18) describes a circle of radius r centered 
A |A I 

at H(f) as pictured in Figure F-3. In terms of the gain factor estimate IH I 


The use of multiple estimates is employed in the fast Fourier transform 
(FFT) type of spectral analysis; more generally the term n^ can be replaced 
by the B^T product, where B is the effective analysis bandwidth and T is 
the total duration of the record. See Equation (F-47). 



and the phase factor estimate «;b, if r < |Hj then the "approximate” 68% 
confidence intervals for the actual gain factor )Hj and phase fa ;tor ^ are 
given by 


iAl - V 


where 


.A - 1 A , , A , 

A<^ = sin (r / 1 H |) 


(F-20) 


(F-21) 


/s, 

r 



Figure F-3. Confidence Interval for Transfer 
Function Estimates 
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The standard error in these results can be estimated by € where 

t . ?/| HI » 

using Equation (F-9) for |H|. For small €, tc obtain the "approximate" 95% 
confidence intervals for gain and phase, the radius r should be replaced by 
2r in Equations (F-20) and {F-21). 

F.4 ESTIMATION OF NOISE -TO -SIGNAL RATIOS USING 
EXTERNAL EXCITATION 

From x(t) = u(t) + m(t) and y(t) = v(t) + n(t), when u(t) and v(t) are due 
to i(t), and m(t) and n(t) are independent of i(t), it follows that 


G. = G. +G. = G. = A. G.., G. = G. + G. = G. = B, G.. 

IX lu im lu lu u ly iv in iv iv ii 


(F-23) 


where A, and B. would be computed transfer functions. The quantities G 

lU IV ^ uu 

and G then become 
vv 


G 

A. 

^ G.. 

2 

= Y G , 

G 

B. 

uu 

lU 

11 

IX xx’ 

VV' 

IV 


^ G.. = G (F-24) 

11 ly yy 


Now G = G + G and G 

XX uu mm yy 


G + G give 
vv nn ® 


mm 


= (i-y^ )( 

\ IX' 


XX 


G = ( 1 -y? )g 
nn V ly ' yy 


(F-25) 


Hence the noise-to-signal ratios 



yff = 


G /G 
nn vv 



(F-26) 
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1 


Conversely, Equation (F-26) proves that 


= l/(l + a), (F-27) 

Note that a = B if and only if = y? . 

^ IX ly 

2 

Values of y. and y. for various values of a as might be used in 

IX IX ® 

Equation (F-22) are shown in Table F-1. 

Table F-1. y^ versus a [Eq. (F-27)] 


a 

y? 

IX 


0.05 

0.952 

0.976 

0. 10 

0.909 

0.953 

0. 15 

0.870 

0. 933 

0.20 

0.833 

0. 913 

0.25 

0.800 

0. 894 


i: . 6 SEl^r iVCISE CASE 

The self-noise case corresponds to i(t) = 0 with u(t) = 0 and v(t) = 0. 
For this case, x(t) = m(t) only and y(t) = n(t) only. Here 


G=G , G=G, G=G 
XX mm yy nn xy mn 


Equation (F-1) now yields 



(F-28) 


(F-29) 


From Equation (F-6) 


H = H - G /G 
xy mn mn mm 


with associated 



and d> . 
^mn 


(F-30) 
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These estimates of y , iH I and d> should be compared to 

2 mn^ I mn| ^mn ^ 

estimates of | ^xy | ^xy whc;n signals are present due to i(t). 

In particular, poor estimates of H v/ill be obtained if y^ is small. 

^ * mn mil 

F. 6 HYPOTHETICAL NO-> NOISE CASE 

In this hypothetical case, it will be assumed that m(t) = 0 and n(t) = 0 
so that x(t) = u(t) and y(t) = v(t). Here, Equation (F-1) yields 


y = y 

' xy uv 

2 

with y^^ = 1 under ideal conditions. From Equation (F-6) 


G I Y' 

uv| / 


G G 
uu vv 


(F-31) 


H = H = G /G 
xy uv uv uu 


(F-32) 


with associated I H 


uv 


and <(> 


uv 


2 

These estimates of |H 1,0 and y =1 will be considered to be 

I uv| ^uv uv 


the desired (unknown) true values of 

F. 7 GENERAL CASE OF CORRELATED NOISE 


It is assumed here that the noise terms m(t) and n(t) are partially 

correlated with each other but not with the signals u(t) and v(t) that are due 

to i(t). This case is covered by Equations (F-1) through (F-5). It is reasonable 

to assume that the two terms in G = G + G satisfy 

xy uv mn ' 


mni 



(F-33) 


Bounds for y^^ of Equation (F-5) will now be developed. 

Two special situations can occuv corresponding to being 

parallel to j or belnj; pe rpendLcuiar to 
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C&8© A# l^mnl to 

For this situation 

I G I = |G + G 1 = 
I xy I I uv mn | 

wi‘h d> = d> as shown in Figure F-4, 

^xy ^uv ® 



G I ^ |G I 
uv| I mn I 


(7-34) 


Figure F-4. G I Parallel to G 
® mn| uv 


It follows that 


G ^ = G ^ 
xy uv 


[iaIg |/|g |]^ = g G [l±y (a)S)^^^l 

I I mnj/ I uv|J uu vv ^ I 


by setting 


mn 


2 .,2 


r * 


L. 

(F-35) 

G , 

> 1 G 

2 

= G G 

(F-36) 

1 nr 

1 uv 

uu vv 


= 

G /G 
nn vv 


(F-37) 
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Now, the coherence function of Equat on (F-5) becomes 




l + a)(l+iff) 


(F-38) 


The maximum reduction in measured coherence will occur by choosing the 
minus sign in the numerator and by letting approach one. Assuming 

<* = this will yield the lowest possible results 

2 


y^y » (l-a)2/(l + a)^ 


(F-39) 


Values of V and y for v^irious values of a as might be used in 
xy xy ® 

Equation (F-22) are shown in Table F-2. 

Table F-2. Minimum y^^ versus a [Eq. (F-39)3 


a 

y 

*xy 

y 

'xy 

0.05 

0.819 

0. 905 

0. 10 

0.669 

0. 818 

0. 15 

0. 546 

0. 739 

0.20 

0.444 

0.667 

0.25 

0. 360 

0.600 


Case B. | perpendicular to 

For this situation 


hxy|' ' |°uv*°mn!' ' |°uv|'- |°m„| 


(F-40) 


with (A = <A ± as shown in Figure F-5. 

xy uv mn ® 
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) 



Figure F-5. j Perpendicular to 

By using Equations (F-3b) and (F-37), it fellows that 

(F-41^ 


(F-42) 

As before, the maximum reduction in measured coherence will occur by 
letting approach one. Assuming a = fi, cnis will yield the results 

« (l-a^)/(l+a)^ (F-43) 

2 

Values of y and y for various values of a as might be used i i 
xy xy ® 

Equation (F-22) are shown in Table F-3. 


I G I ^ = G G r 1 - y^ 
i xyl uu vv L ni 


Hence the coherence function y^^ of Equation (F-5) becomes here 
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Table F-3. Expected versus a f Eq. (F-43)] 

y' y 

g xy 'xy 

0.05 0.905 0. 951 

0. 10 0.818 0.905 

0. 15 0.739 0. 860 

0.20 0.667 0. 817 

0.25 0.600 0. 775 


F. 8 STANDARD ERROR IN TRANSFER FUNCTION 
ESTIMATES 

Return now to Equation (F-22) which gives the standard error € in 

A 

transfer function estimates H that are obtained by using the external excita- 
tion signal as given by Equation (F-13). Maximum values of C will occur 

2 2 2 
when y and y. take on their minimum values. This will occur for y 
xy IX 2 

from Equation (F-39) and for y. from Equation (F-27), Substitution of these 

XX 

results into Equation (F-22) gives 


€ < 


L 

(1-g) 


g (l+g) 
"d 




1/2 


(F-44) 


indicating how € varies as a tui^ction of n^ and a. 

To illustrate use of Equation (F-44), suppose a = 0. 10, corresponding 
to y? « 0, 909 and y^ ^ 0, 669 from Tables F- 1 and F-2. Now 

' XX xy 


€ 


1 0-74 


Hence n, = 25 gives £ S 14. 7% and n , = 100 gives £ S 7.4%. In terms of £, 
a a . 

there will be about 68% confidence that measured values of |H | will be within 
* ^% of the desired unknown true values | H|, and there will be about 95% 
confidsnce that measured values of |H| will be within * 2£% of the desired |H|. 
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Note that for small values of a. Equation (F-44) can be approximated 


by 



(F-45) 


F. 9 RELATION OF NUMBER OF AVERAGES TO 
NOISE^TO-SIGNAL RATIOS 

The number of n^ of independent averages can be related to the total 
available record length formula 


T^ ^ = n ,T 

total d 


(F-46) 


where the subrecord length T = l/B^^ the quantity being the bandwidth 
resolution. Hence 


n, T^ ^ ,/T = B T^ ^ , 
d total e total 


(F-47) 


Multiplying numerator and denominator of the expression for the noise -to- 
signal ratio a in Equation (F-26) by B yields 


a = G B 
mm e 




2G B /u 
mm e/ 


(F-48) 


2 . 


where u is the mean square value of u(t) and u^ is the amp ide of u(t) if 

it is a sinusoid. 

Substitution of Equation (F-47) into Eqixation (F-48) shows 


a = G 


mm 


"d/" 


^ , = 2G n,/u^T. . , 
total mm d# o total 


(F-49) 


Thus, a increases as n. increases and a decreases as T^ ^ - increases. 

d total 

Also, Equation (E-45) becomes now 

11/2 


€ S 2 


n 1*'^ r iz 1 

total J = ^ 


1/2 


(F-50) 


where would be measured by Equation (F-25). 
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